









Volume 19, No. 4 





ABTOMATHKA April 1958 | 

TEAEMEXAHKMKA | 

Sk. ne OF MICHIGAN | 
MAY 4 - 1959 


".= 


SOVIET INSTRUMENTATION AND 
CONTROL TRANSLATION SERIES 


Automation 
an 
Remote Control 


(The Soviet Journal Avtomatika i Telemekhanika in English Translation) 


@ This translation of a Soviet journal on automatic control is 
published as a service to American science and industry. It is 
sponsored by the Instrument Society of America under a grant 
in aid from the National Science Foundation, continuing a pro- 
gram initiated by the Massachusetts Institute of Technology. 

















Instrument Society 
of America 
Executive Board 


Henry C. Frost 
President 


Robert J. Jeffries 

Past President 

John Johnston, Jr. 
President-Elect-Secretary 
Glen G. Gallagher 

Dept. Vice President 
Thomas C. Wherry 
Dept. Vice President 
Philip A. Sprague 

Dept. Vice President 
Ralph H. Tripp 

Dept. Vice President 
Howard W. Hudson 
Treasurer 

Willard A. Kates 
Executive Assistant—Districts 
Benjamin W. Thomas 
Executive Assistant—Conferences 
Carl W. Gram, Jr. 

Dist. I Vice President 
Charles A. Kohr 

Dist. II Vice President 

J. Thomas Elder 

Dist. III Vice President 


George L. Kellner 
Dist. IV Vice President 


Gordon D. C e. 
Dist. V Vice pe a 
Glenn F. Brockett 

Dist. VI Vice President 
John F. Draffen 

Dist. VII Vice President 
John A. See 

Dist. VIII Vice President 


Adelbert Carpenter 
Dist. 1X Vice President 


Joseph R. rs 
Dist. X Vice President 


Headquarters Office 
William H. Kushnick 
Executive Director 

Charles W. Covey 

Editor, ISA Journal 

George A. Hall, Jr. 

Assistant Editor, ISA Journal 
Herbert S. Kindler 

Director, Tech. & Educ. Services 
Ralph M. Stotsenbu 

Director, P seoters sy 
William F. Minnick, Jr. 
Promotion Manager 





SOVIET INSTRUMENTATION AND 
CONTROL TRANSLATION SERIES 


ISA Publications Committee 


Nathan Cohn, Chairman 

Jere E. Brophy Richard W. Jones John E. Read 

Enoch J. Durbin George A. Larsen Joshua Stern 

George R. Feeley Thomas G. MacAnespie Frank S. Swaney 
Richard A. Terry 


Translations Advisory Board 
of the Publications Committee 


Jere E. Brophy, Chairman 
T. J. Higgins S. G. Eskin G. Werbizky 


@ This translation of the Soviet Journal Avtomatika i 
Telemekhanika is published and distributed at nominal 
subscription rates under a grant in aid to the Instrument 
Society of America from the National Science Foundation. 
This translated journal, and others in the Series (see back 
cover), will enable American scientists and engineers to 
be informed of work in the fields of instrumentation, 
measurement techniques and automatic control reported 
in the Soviet Union. 


The original Russian articles are translated by competent 
technical personnel. The translations are on a cover-to- 
cover basis, permitting readers to appraise for themselves 
the scope, status and importance of the Soviet work. 


Publication of Avtomatika i Telemekhanika in English 
translation started under the present auspices in April 
1958 with Russian Vol. 18, No. 1 of January 1957. Trans- 
lation of Vol. 18 has now been completed. The twelve issues 
of Vol. 19 will be published in English translation by 
mid-1959. 


All views expressed in the translated material are intended 
to be those of the original authors, and not those of the 
translators, nor the Instrument Society of America. 


Readers are invited to submit communications on the qual- 
ity of the translations and the content of the articles to 
ISA headquarters. Pertinent correspondence will be pub- 
lished in the “Letters” section of the ISA Journal. Space 
will also be made available in the ISA Journal for such 
replies as may be received from Russian authors to com- 
ments or questions by American readers. 


Subscription Prices: 
Per year (12 issues), starting with Vol. 19, No. 1 


General: United States and Canada. . .. . $380.00 
SG at eeke ar we 6 se 
Libraries of non-profit academic institutions: 
United States and Canada. . .. . $15.00 
tg De sg = gw, ae oe. aoe 
Single issues to everyone,each. . ... . . $ 6.00 


See back cover for combined subscription to entire Series. 


Subscriptions and requests for information on back issues 
should be addressed to the: 

Instrument Society of America 

313 Sixth Avenue, Pittsburgh 22, Penna. 


Translated and printed by Consultants Bureau, Inc. 








Dr. 
Dr. 
Dr. 
Dr. 
Dr. 
Dr. 


Sci. 
Sci. 
Sci. 
Sci. 
Sci. 


Sci 


See following page for Table of Contents. 





Volume XIX No. 4—April 1958 





Automation and Remote Control 


The Soviet Journal Avtomatika i Telemekhanika 
in English Translation 


Avtomatika i Telemekhanika is a Publication of the Academy of Sciences of the USSR 


EDITORIAL BOARD 
as Listed in the Original Soviet Journal 


Corr. Mem. Acad. Sci. USSR V. A. Trapeznikov, Editor in Chief 
Dr. Phys. Math. Sci. A. M. Letov, Assoc. Editor 

Academician M. P. Kostenko 
Academician V. S. Kulebakin 
Corr. Mem. Acad. Sci. USSR B. N. Petrov 
Tech. 
Tech. 
Tech. 
Tech. 
Tech. 
Tech. ‘ 
Cand. Tech. Sci. V. V. Karibskii 
Cand. Tech. Sci. G. M. Ulanov, Corresp. Secretary 
Eng. S. P. Krasivskii 
Eng. L. A. Charikhov 


M. A. Aizerman 
V. A. Ilin 

V. V. Solodovnikov 
B. 8. Sotskov 

Ia. Z. Tsypkin 

N. N. Shumilovskii 


Copyright by Instrument Society of America 1959 

















English Translation Printed February, 1959 





AUTOMATION AND REMOTE CONTROL 





























































Bibliography 
List of Existing Literature on Magnetic Amplifiers and Contactless Magnetic Components . . 


Volume 19, Number 4 April 1958 
RUSS. 
PAGE PAGE 
Graphico-Analytic Method for Determining Relay System Characteristics, L. P, Kuz'min.. 277 285 
Concerning a Method for Analyzing Sampled-Data Systems, Fan Chun-Wui,.......... 288 296 
On Improving the Transient Response of Correcting Links with Variable Parameters. E. K. 
DPCM EEE Cab a Gemeccccoceccceceecoeeeseeesecsooeeocte 299 306 
Statistical Investigation of Nonstationary Processes in Linear Systems by Means of Inverse 
Simulating Devices, EEE a ar e- erer aea 305 312 
The Influence of Fluctuations on the Operation of an Automatic Range-Finder. I. N. 
Amiantov and VY, I. Tikhonov eeeeeeeeeeeeeeeeeeeeveeeeeeeeeeeeeeeee 318 325 
Determination of System Parameters from Experimental Frequency Characteristics, A. A. 
i i . occ c ease ede eee ceeeecceeeeseeeees 327 334 
Optimal Frequency Deviation in One-Channel Telemetering System, Iu, 1, Chugin...... 339 346 
A Selective Low-Frequency RC-Amplifier as a Control System Element. Iu, G, Kochinev.. 349 355 
Phase Detector for Multiple Frequencies, R. Ia. Berkman. ...ccecceccccceseseces 355 360 
Concerning Some Properties of Ferromagnetic Clutches, P, N. Kopai-Gora .....+e+++++ 361 366 
ty Additions to the Table of Optimal Characteristics Given by Solodovnikov and Matveev. 
: Dis Mee IE Be So cece cer ececccecoe ee eesecoecconcoreececeses 370 376 

































/2~/492 ¢¢ 


2,0 


GRAPHICO-ANALYTIC METHOD FOR DETERMINING RELAY SYSTEM 
CHARACTERISTICS 


L. P. Kuz‘min 


(Moscow) 


A graphico-analytic method is presented for the determination of 
relay system characteristics J,;(w) and },(w), necessary in the investiga- 
tion of periodic regimens in relay systems when the frequency method is 
used, There will be given an example of the determination of the 
characteristics J,(w) and J,(), using the actual frequency characteristic 
of the linear portion of the system. 


In the investigation of periodic regimens inrelay designs, the schematic diagram is usually given in the 
form shown in Fig. 1. All the linear (linearized) links are treated as one link, characterized by the transfer 
function Wjp (p). The external controlling function u(t) and the external stimulus f(t) are brought to the input 
of the relay element, In this case, the controlling signal, x(t) at the relay input is made up of the controlling 
function u(t), f(t) and the feedback signal 2{t). As a result of the impression of the controlling signal x(t), a 
signal y(t) is produced at the relay output, this signal taking the form of a sequence of rectangular pulses of 
constant amplitude, Acting on the linear link, these signals give rise to the appearance of the output signal z(t) 
of the closed-loop relay system. 


If the system is in a periodic regimen, the quan- 
tities x(t), y(t) and z(t) vary periodically with a definite 
frequency, The investigation of the periodic regimen 
in relay systems boils down to the determination of the 
periodic solution of the equation 


Z(t) = Wip (p) y (2), * (1) 


satisfying the conditions that the given periodic regimen a 
exist [1]. anes 





y(t) , so) 2() 
































Fig. 1. Schematic diagram of a relay 


In the majority of cases, the direct determination of the periodic solution of Equation (1) is a quite difficult 
problem, Therefore, it is generally solved by means of various special methods [2-5]. 


The investigation of periodic regimens in relay systems by the frequency method [2] is carried out by means 
of the so-called characteristic of the relay system 


J (wo) = —27 (2) = jz). (2) 


If the relay element possesses an insensitive zone, then the periodic regimen is characterized not only by 
the frequency w, but also by the relative time during which the relay contact is in the closed state, y = t,,/T, 


* Here and in the sequel, the sign ~ over the symbol for a quantity indicates that the quantity varies periodically. 








te sx 


where T is the half-period of the given periodic regimen. 
In this case, the investigation of the relay system is carried out by means of two characteristics: 


i</n\ ~/% 7“/ gs , r\ 
—=* (S)—#(F)s 4 =—52 (1 E)—*y eS), (3) 
which are functions of: two parameters: the frequency w and the relative time the contact is in the closed state, 
Y- 


The characteristics (2) and (3) of the relay system are closely related to the transfer function of the linear 
portion, Wip (p), and, in those cases where the poles of Wj)(p) are known, may be calculated from this transfer 
function, Benseven, J(w), Jy ) and Jy(w) may be cocmmied from the transient (timewise) and from the amplitude~- 
phase characteristics of the linear portion [1]. 


The latter method is used the most frequently since, in practice, it has no limitations from the point of view 
of structure of the system, parameters, stability of the linear portion, etc. It provides the possibility of calculation 
of the characteristics (2) and (3) from empirically determined amplitude-phase characteristics of the systems 
being investigated, 


The relay system characteristics are connected to the amplitude~phase characteristic of the linear portion 
by the following expressions, given in the form of series: 


J, (o) = —_ > [ {1 —cos(2m—1) yx]u{(2m—1)<0]—sin(2m—1 )yxo[ (2m—1 jo]} + 


m=1 








$+ {ASPX wu f(2m — 1) 0] + SET yp ((2m — 1) ol], 
Jy (@) = mi > [ {(4—cos(2m—t )y]uf(2m—1 eo) + sin(2m—1)yxv[(2m—1)o]}— (4) 
m=1 








—j (eon ee u [(2m — 1) @] — =e {YF 5 (2m — 1) a]}], 


where u(w) and v(w) are, respectively, the real and imaginary parts of the amplitude-phase characteristic. 


If there is no intensitive zone (y = 1), then the characteristics J;(w) and Jy(w) coincide, and are determined 
from the following expression: 


ae v cA) 
J(o)=—2 J [ut(@m— 1)@] +7 a) il. (5) 


m=1 





Analysis of Expressions (4) and (5) shows that the determination of the relay system characteristic and, in 
the former case, of the characteristics Jw ) and J y(w), entails a great deal of computational work. In fact, to 
determine J,(w) or by (w) from Formula (4) it is paeiaaey. for each fixed value of the frequency, to compute mp 
terms of each series (here, m is the number of terms taken, and p is the number of fixed values of the parameter 
Y). There is thus the constant necessity to compute mpq terms in Expression (4), where q is the minimum quan- 
tity of fixed frequencies necessary for the construction of the characteristics Jy(w) and Jy(w). 


The magnitude of the computing job essentially complicates the use of the frequency method in investiga- 
ting relay systems, This is particularly the case when a choice of parameters for the periodic regimen is being 
made, with concomitant variations of the structure of the design (introduction of feedback, etc,), since here the 


expressions J,(w) and Jy(w) are repeated more than once. 


The problem of investigation is essentially simplified if a graphico-analytic method is used for determi- 
ning the characteristic, In essence, the graphico-analytic method of computing the characteristics J,(w) and 
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Jy(w) consists of picturing their components, corresponding to the values m = 1, 2,..., 00 and y = 1, in the form 
of vectors, for each fixed value of frequency, determining the hodographs of these vectors when y varies within 
the limits 0 <y <1, and summing the vectors for equal values of the parameter y. The resulting vectors deter- 


mine the points of the characteristics J;(w) and Jy(w) forthe given frequencies and for all values of y in the 
interval O< y <1, 


As will be shown below, the hodographs of these vectors are determined by an uncomplicated geometric 
construction, which also gives the possibility of facilitating the determination of the relay system characteristic, 


As a basis for our method, we transform Equation (4). Considering that, in practice, the determination of 
J;(w) and Jy(u) is usually restricted to n terms of the series, we present Expressions (4) in the following form: 


J, («) = Ju, (0) + Ja, (30) + «++ + Sa, ((2n—1) a0), 
Jy (00) = Ty, (0) + Syg (30) +++ + Jy, (22 —1) 0). 


Here, the complex terms J; (w), 31,(3w), . oes ly Phe J -(3W)- + corresponding to the values m = 1, 2,..., 
n, are determined by hetenteg te corresponding vA ues of'm in the expression 


Tay, ((2m — 1)«) = on = [fi — cos (2m — 1) yx} u((2m — 1) w) — 





— sin (2m — 1) yxv (2m — 1)«)} + j {GRO u((2m— 1)0)+ 


— ee 


+ — v ((2m — 1) «)}], 


Tem ((2m — 1) w) = * lu — cos (2m — 1) yx} u((2m — 1) 0) + 
+ sin (2m — 4) ynv ((2m — 4) o)} j {=O oy ((2m— 1) a) — 


2m 








__ 1— cos (2m — 1) yx 


mT? ((2m—1) o)}]. 


The expressions in (7) correspond, in the complex plane, to vectors, the magnitude and direction of which 
are determined by the frequency w and the relative time of contact closure y, but for fixed values of the fre- 
quency, w = w,, depend only on the parameter y. With continuous variation of y within the limits 0<y<1, the 
ends of each of these vectors describe hodographs, The analytic expression for the hodographs of the vectors 
Jim andJ,,, may be determined immediately from Equations (7) by substituting in them the fixed value of 
frequency, w = Wye 





Transforming Formulas (7), we present the analytic expression for the hodographs J;,,, and Jy, in a clearer 
form, 


First, we consider the hodographs of the vectors J;,,, and Tym for m = 1, In this case, substituting in Expres- 
sion (7) m = 1 and w = w, and carrying out a simple transformation, we obtain 


S13, (Y%) @x) = > tu (cox) + 7 (wx)] — [u (@x) + Jv (@x)] (cosyx — 7 sin y)}, 


®) 
Jy, (Xs 4) = =P {fu (ox) + 70 (on)] — [w (coe) + 70 (@)] (cos yx + j sin yn). 
Letting 
£ (om) = =P (w(x) + jo(ox)) = —2W (jon), e 
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Fig. 2. Hodographs of the vectorsJ;__, [(2 m —1)w,] andJy,, [((2m —1),] 

for m = 1, 
and also taking into account the Euler equations, we can write the equations for the hodographs of (7 » W,) and 
Jy ’ wy) in the following form: 


Ja, (Yr @e) = 2 (@e) — 2 (Ox) e-P™, Ty, (Y, Oe) = 2 (@~) — 2 (x) CM. (8") 


An equation of the type (8) or (8") is the equation of a circle of radius z(w,) whose center has been dis- 
placed to the point z(w;,). Bearing in mind that in Equations (8") the argument varies from 0 to #, and that the 
displacement from its center to the origin of coordinates is equal to its radius, we see that the hodographs of 
Jig(7 » Mand Jy, (7, W) are semicircles, cutting the origin of coordinates when y = 0 (Fig. 2). The sum of the 


hodographs of Iu, CY» Wy) and Jy (7, wy) Isa full circle whose diameter is 4k)/T, the modulus of the ampli- 
tude-phase characteristic of the linear portion of the system. In fact, we find from Equation (8) that, for y = 1, 


a ee 4k 
Jay (Ys x) = Sy, CY) On) = —® W (jou) = = W, (con) ee, 
where W,(w,) is the modulus of the amplitude-phase characteristic when w = wy. 


By analogous transformation, we can put the equations for the hodographs Jip (1s Wy) and Jy ml? » Wy)» for 
m > 1, in the form 


J smn, (Ys (2m — 1) oy) = 2, — =n ei(2m—1)rn __ = e—i(am—1)re 
Tn (Ys (2m — 1) top) = Ziq — ABFA ftem—ayrn Fim = Fim sammy, (10) 








where 





2k 2m — 1 
2105 =u ((2m — 1) wy) +7 hs! |, 


2k, u ((2m — 1) ,) 


| qa 
tm = | wae + jo((2m—1)oy)|. ’ 





Equations (10) represent ellipses with centers displaced to the point z4,,, (Fig. 3), determined from Expres- 
sions (11), 


It is easy to see from Equations (10) that the hodographs of vectors Jy,(7 , w,) and jy mm » Wy) coincide, 


differing only in the direction of rotation of the vectors J 4,, and} ,-* 


Jim}, 


on 
=? V/sw) 
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Fig. 3. Hodographs of the vectors 3,,,{(2m —1) uJ Fig. 4. Construction of the hodographs for the vectors 
and Jy,,{(2m —1) w,] for m = 2, J imf(2m —1) w,] and Jy mi(2m —1)w,] for m = 2. 


For m = 1, the equations for the ellipses, (10), 
degenerate into those for the circles (8"), Considering 
circles as particular cases of ellipses whose semi-axes 
are equal, we are led to the conclusion that the hodo- 
graphs of the vectors J;,, and Jyn are ellipses, the 
magnitude and position of which are determined, for 
fixed frequency, by the magnitude of the real and imagi- 
nary parts of the amplitude-phase characteristic of the 
linear portion of the system, 


It follows from this that, in the determination of 
Fig. 5. Determination of the vectors Wij (2m —1) JxW) and Jy(), it is not necessary to compute, for every 
wk] from the amplitude-phase characteristic of the value of the frequency, the magnitude of the components 
linear portion, J1,, 2nd Jy, for different values of the parameter y. It 
is sufficient just to find the components of J;,,, and Jy 
for y = 1, and all the remaining ones, corresponding to other values of y, may be found from the hodographs 
determined by Expressions (10) and (8"). 


jVlw) 











The construction of the hodographs of the vectors Jj, and J,,,, and the determination of the points of these 
hodographs corresponding to different values of y, is carried out according to Equations (10), However, it is 


*When y varies between the limits 0 < y = 1, the vectorJ,,, rotates clockwise, and the vector Jy, rotates coun- 
terclockwise (Fig. 3). 
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simpler to carry out these operations graphically. 


In analyzing Expressions (7) for m> 1, and in comparing them with the corresponding expressions for m = 1, 
we note that, for m > 1, there is a decrease in the imaginary parts of Jip, and Jy, and also an increase in the 
arguments of trigonometric functions sin and cos by a factor of (2m—1). The first condition entails a deformation 
of the hodograph in amplitude, and the second entails an increase in the argument of the vector Jy, dg, with a 
change of the parameter y in comparison with the increase of the argument of the vector J; 1%, ), again by a 
factor of (2m —1). 


Bearing these in mind, and also taking into account that, for m = 1, the hodograph of the vectors J,,,, and 
Tym 18 4 circle, we establish the following sequence for constructing the hodographs of these vectors for m > 1: 


a) on the complex plane the amplitude-phase characteristic corresponding to the frequency (2m —1) w is 
constructed, magnified by the factor 4k, / TT; 


b) with this vector as a diameter, an auxiliary circle is constructed which is then divided into 2/k parts, 
where k is the chosen interval of variation of the parameter y (the origin of coordinates is taken as the starting 


point); 


c) for the determination of the vectors Ji,, and Jy, corresponding to various values of y, we choose 
various values of this parameter and, using the markings on the auxiliary circle, we count off from the origin of 
coordinates (2m —1) y /k intervals, counting clockwise for the vector J;,,, and counterclockwise forJ,_. Decreas- 
ing by a factor of (2m —1) the ordinates of the points thus obtained, we get the points of the hodograph for Jin 
Gy,,) corresponding to the selected values of y. 


Figure 4 gives the construction of the vectors },,,, and Jy, for m=2 and w = 3w, = w; for two values; 
: = 0.1 and y = 0.5. Vectors OA and OB correspond to the vector hh, (3w;,), and OC and OD to the vector 
Jy (3u,). For y = 1, vectors Ji, (Sw) and yy, (3u;,) coincide (vector OZ). The dotted lines show the hodograph 


fof the vectors Ty, (3wy) and By (3w,): 


In determining the amplitude-phase vector necessary for starting the construction of the hodographs, one 
can make immediate use of the amplitude-phase characteristic, reading off the points on it corresponding to the 
frequencies w,, 3wy, 5wy,...(2m—1) w,. The vectors from the origin of coordinates to these points are pro- 
portional to J; (wy), Jg(3w,), Jg(5uy,), ..- Imm [(2m —1) wy, ] (Fig. 5). 


It is more convenient to carry out the determination of the starting vectors from the family of curves for 
the real, uf{(2m —1)w], and imaginary, v[(2m —1) w], parts of the amplitude-phase characteristic (Fig. 6), from 
which the range of variation of m for all computed frequencies may also be determined. The family of curves 
u{(2m —1)w] and v[{(2m —1)w] are constructed from the real, u(w), and imaginary, v(w), parts of the amplitude- 
phase characteristic by recomputing the arguments. By setting the values-of the frequencies we then, according 
to the ordinates of the corresponding curves (u(w)and v(w), u(3w) and v(3w), etc.), find the vectors of the ampli- 
tude-phase characteristic, Multiplying these by the factor 4k,/m , we obtain the initial vectors J, (wy), Jz (3), 
-+++» Bmf(2m —1)u,]. After this, we determine the range of variability of m. According to Fig. 6, one can 
limit ones consideration to the components of J,(w) and J(w) corresponding to m = 1 (i.e., the range of variabi- 
lity is m = 1). For frequencies wy < w < ws one can limit oneself to m = 2, for frequencies w, < w < Wy to 
m = 3, etc. 


The construction of the curves u{(2m —1) w] and v{(2m —1)w] is simplified if the logarithmic scale is chosen 
for the frequency [1]. Bearing in mind that one octave corresponds to a doubling of the frequency, these curves 
are simply constructed by translating the curves of u(w) and v(w) to the left 1.6 octaves for m = 2, 2.3 octaves 
for m = 3, etc. (Fig. 7). 


Taking all of the above into consideration, we can set up the sequence of steps for the determination of 
the characteristics J;(w) and }y (w): 


1) the amplitude-phase characteristic of the linear portion of the system is computed, and the curves u(w) 
and v(w) are constructed; 


2) from the curves u(w) and v(w) are constructed the families of curves u{(2m —1)w] and v[(2m —1)w] 
for m = 2,3,...,0; 








u(( 2m-Nw) 














Fig. 6, Family of curves u{(2m —1) w] and v{(2m —1)w]. 
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Fig. 7. Construction of the curves u{(2m —1)w] and v{(2m —1) w] for 
m = 2,3... from the curves u(w) and v(w). 
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Fig. 9. Hodographs of the vectors Jy,,. [(2m —1) w,] and Tym [(2m —1) w,) and 
the graphic determination of the vectors (7 ,w,) and Jy (y ,w,) for wy = 20 
radians/sec, and m = 2, 


3) the range of variation of m is determined for all computed intervals of the frequency w; 


4) for the chosen computed values of the frequency, w = wy, the initial vectors J,,[(2m —1)w,) = 
JT %p[(2m —1) Wy] and their hodographs are constructed; 


5) we add the vectors J;,.[(2m —1)w,] and Jy mi(2m —1) w,] which correspond: to m = 1,2,...,n and to the 
same value of the parameter y, thus obtaining the vectors J;(w,) and J (w,) for the given frequency and for 
various values of y in the range 0< y= 1; 


6) similarly, we determine the vectors J y(w) and Jy(w) for all other computed v@lues of the frequency, after 
which we construct from them the characteristics J,(w) and Jy(w) [or their imaginary parts, ImJq(w) and ImJ,(w)], 
necessary for the investigation of periodic regimens in relay systems, 


Thus, using the graphico-analytic method for determining the characteristics of relay systems involves the 
computation of just one amplitude-phase characteristic, The remaining operations necessary for the determination 
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Fig. 10. Hodographs of the vectors J;(y, ws) and 
3, (7 Wy) for ws = 40 radians/sec, and m = 1, 








Fig. 12. Hodographs of the vectors J;(y, ws) and 
Jy (7 ,Wg) for ws = 80 radians/scc, and m = 1, 








Fig. 14, The curves ImJ, (w). 
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Fig. 11. Hodographs of the vectors J; (y , wg) and 
ie ,W,4) for ws = 60 radians/sec, and m = 1, 
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Fig. 13, The curves ImJ; (w). 


of the characteristics J,(w) and Jy(w) are carried out 
graphically, which significantly simplifies the problem of 
determining these characteristics and, consequently, the 
entire problem of investigating periodic regimens in relay 
systems. 


As an example, we compute the characteristics 
J,x(w) and Jy(w) for the relay system for which the families 
of curves u[{(2m —1)w] and v[(2m —1) w] are shown in 
Fig. 7. 


We determine the range of variation of mfor various 
values of the frequency: m = 1 for w = 30 radians/sec.; 
m = 2 for 20 radians/sec. sw < 30 radians/sec.; m = 3 
for 15 radians/sec. sw < 20 radians/sec,; m = 4 for 
w< 15 radians / sec, 


We set the values of the frequency: wy, = 10 
radians/sec, for m = 4, w, = 20 radians Aec. for m = 2; 
Ws = 40 radians/sec., wg = 60 radians/sec., ws = 80 
radians/sec, for m = 1. 


We set the values of the parameter y: y, = 0.2, vy, = 0.4, 73 = 0.6, yg = 0.8, yg = 1.0. 


Using the method described above, we determine the starting vectors Jy,,[(2m —1) w,] and construct the 


hodographs (Figs. 8-12). 














Summing the vectors corresponding to one value of the parameter y, we find Jx(y , Wj) and J, (7, Wy) for 
the frequencies wy and wy (Figs, 8a and 9a). For the frequencies ws, wg and ws, for which m = 1, the vectors 
3x(Y, w,) and Jy (7, W,) are determined directly from the hodographs, shown in Figs. 10-12, On Figs. 8a, 9a, 
10, 11 and 12, the solid lines represent the vectors J;(w,), the dotted lines represent ¥(w,). 


The resulting vectors thus obtained determine the points of the characteristics J,(w) and Jy(w) for various 
values of the parameter y. On Figs, 13 and 14 are shown the imaginary parts ImJ,(w) and IMJ(w) of these 
characteristics for values ReJy (w) < 0 and ReJ,(w) < 0, which are necessary for the direct determination of the 
existence of the periodic regimens and their parameters, 
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CONCERNING A METHOD FOR ANALYZING SAMPLED-DATA SYSTEMS 


Fan Chun-Wui 


(Moscow) 


Based on the discrete Laplace transform, a well- 
known method [1, 2] for analyzing continuous linear 
systems is generalized to sampled-data systems and ap- 
plied to the qualitative prediction of their forced motion. 


It is well known that the output coordinates of a linear system are completely determined by the distribu- 
tion of the zeroes and poles of the closed-loop system's transfer function, the initial conditions and the external 
forcing function. The applicability of such an approach to continuous linear systems was asserted by S. P. 
Strelkov [1, 2]. In articles [3-5] attempts were described to investigate sampled-data systems in an analogous 
manner. However, the case of an arbitrary forcing function was not investigated there. This paper considers the 
accuracy of a sampled-data system to which a discrete forcing function of a very general form is applied. 


Preliminary Remarks 





We consider the sampled-data control, or regulating, system, the schematic for which is shown in Fig. 1, 
where f [n], E{n] and x{n, « ] are some step functions, 


The transfer function for the closed-loop sampled-data 









































ln} Ela] | Pe system has the following form [6]: 
° X34, 
z[n,e] ®;;(q, e) = ae = 
aa hind bot bie +... +b,e" — Ay (4, ®) (1) 
ay + aye? +... + a,c G" (q) 
(j =I, I) 
Fig. 1, 


where xia, € ) and F®(q) represent the discrete transformations of the step functions x{n, € ] and f [n]. 


In Equation (1), m <2, and the quantity ¢€, which can vary within the intervals from 0 to y and fromy to 
1(y is the “porosity"), enter into the coefficients by, by,..., bm. 


We assume that the discrete forcing function f[n] can be represented in the form 
8 
f{n)} = >) ce (ne (2) 


k=1 


where c,[n] is a polynomial in n and p, is a real or complex number. 








[=m 








We write the expression for the output function, x(a. €), in the general form: 





; H;(q, 8) <« H3(q, ©) F;(q) 
ro = Bat Bh. doi =. j LU = 
X3@. =o * M=—@ F(@) 
. ° 3 
H; (q, €) Fy (¢) ” 





(e% — ee — ce)... (eT —eAl) (eF—e™Ye® — et), . , (oF — oP 8) 


Here )4, Ag, ..-, Ay aNd py, Pgs++, Ps are the poles of the transfer function of the system itself and of 
the transformation F*(q) of the forcing function, Initially, for simplicity in the argument, we shall assume that 
all the poles are distinct, Then, by the usual formulae, we find the step function at the output in the following 
form: 


I * . 8 . . 
sin, el = GO) FO) Py G (ey) eRe) 








i=—1 
(4) 
= DS} Siei” + >) Reerk” = xz[n] + xp [n]. 
i k 


Here x; may be called the "inherent" and xp the "forced" components, 


It is clear from Equation (4) that both these components depend on the distribution of the zeroes and poles 
of the transfer functions of the system itself and of the function F*(q). Obviously, the inherent motion of the 
system always entails a necessary error of reproduction, For stability of the system, the real parts of all A; cer- 
tainly must be negative. Moreover, it is clear from Equation (4) that the quantities Ry, are not always proportional 
or equal to the quantities cy of Equation (2). This means that the forced motion of the system is different from 
the controlling signal f[n], i.e., it is distorted. Thus, for accurate transmission of the signals, it is necessary to 
decrease these two components of distortion. 


We shall consider the conditions for a distortion-free forcing part, figuring the inherent motion as having 
a certain error, which will be investigated later, 


Conditions for Distortion-Free Forced Motion 





Initially we shall assume that the inherent motions of the system die out after a certain time, In this case, 
the accuracy of signal reproduction is completely determined by the system's forced motion, We consider the 
following cases, 


1, The poles of the system's transfer function are simple, different from its zeroes, and none of them coin- 
cides with a pole of F®(q). In this case, it is clear from Equation (4) that the sequence of polynomials at the 
discrete moments Rj{n] will equal the sequence of polynomials c,{n], Hence, the conditions for distortion~free 
forced motion will be 


H; (ep, 
—_ = const = C. (5) 
G (p,) 
Condition (5) is not sufficient if F*(q) has multiple poles, For example, if F*(q) has two multiple roots, 
Pig, 1.€., 
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—e Fog) | 
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then f [n] takes the form 


f{n}] =D? [F’ (q)] = 4 (n— 1) "1,2 + 


Fy (P12) 
F 

2 (0,9) ; (6) 
F,(@%) 2 


a F(a) 7 
ef™—1)1,9 4 P) (ek — eP1,.2)* kk, (e, ) ro 
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det Fy (@) qI=°1,9 








forn2= 1, 


The transform of the output function xj (q, €) will be 


H;(q, ©) Fi (9) 
G(q)  Fy(q)(e*— e122)? © 





X3(q, ¢)= (1) 


In this case, the forced motion at the system output will be 


_ Aj (1,9, ©) Fy (P1,2) 
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rl G (P19) F(x) 
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forn = 1, 


It is clear from (6) that, in this case, the conditions for distortion-free forced motion will be 


Hy (x, *) _ Hj (01%, ©) 
G" (ey) G (P;.2) 


ad Fi@ Hj) 
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= const =C, (9) 
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I=?1.9 
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lf Ff (q)/Fy (q) = const, then (10) reduces to an equation of the form 


d H5(a,®) 


——— andy = 0. (11) 
det G (q) 


1.2 





If the poles of F®(q) have multiplicity m(m > 1), then the conditions for distortion-free motion, in ad- 
dition to Equation (5), will be the following equations; 


am™—-1_ F; (9) 
de%m—1 FY (9) 


a™-1 Hi; (q, ©) F; (9) 
detm—1 G"(q) Fy (9) 
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+m (cont'd) 
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deIm—2 G (q) F, (q) Py sm de4m—2 Fy (9) e1+m 
s slauple dicie sethteuk. acanatecshaseee sneeseiieiel ante (12) 
a H5(q, &) Fy (9) d FQ) 
det Gq) Fs (9) |e, 2m de? F)(q) |ey am 
In particular, if Ff (q)/Ff (q) = const, then the System (12) takes the form: 
a™-1 Hi; (q, ®) P 
deim—1 Gq) lem 
am™—2 Hj (a, ©) 0 
detm—2 G°(q) |y.sm (13) 
ad H3(q, ©) my | 
de® ag (q) elm bs : 





Conditions analogous to (13) were obtained in [8], 


2. Poles of F*(q) coincide with poles of the system's transfer function, for example, py = Ay. For simpli- 
city of exposition, we shall assume that the poles of F*(q) and of the system's transfer function are simple and 
different from their zeroes, In this case, the forcing function f[n] may be put in the form 


{{n} = + As F*(q)e™ dq =} Fi) _ yn (14) 
om ch ra ee 


and the forced motion at the system output will be 
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G" (p1) Fy (e1) 
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Opn 


It is clear from a comparison of Equations (14) and (15) that for such systems it is not possible to have 
accurate reproduction of the forcing function f[n]. Therefore, we add to the conditions for absence of distortion 
the condition that there be no resonance, i e., we impose the requirement that the poles of F®(q) do not coincide 
with the poles of the system's transfer function, 


3. The poles of F*(q) coincide with the zeroes of the system's transfer function for errors, The transfer 
function for closed-loop sampled-data systems for errors will be 


+. Bin(@) _ D* (a) 
90 (1) = Fay = a) a6) 





The step function for the errors, E{n], will be 
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(17) 





It is clear from (17) that for accurate reproduction of the effective signals, the system must be such that 
the poles of F*(q) coincide with the zeroes of the transfer function of the system for errors, i.e., that 


D (px) = 0. (18) 
In particular, if F®(q) = Ff @/ cet —1)7, then the system's transfer function for errors must be 


, Dy (q) (e* — 1)* 
®.. (9) = G*(q) ? (19) 





Inherent Motion 





In order to decrease the error in reproduction, it is necessary to shorten the interval of time in which, after 
initiation of the process, there is an appreciable inherent motion. According to (4), the inherent motions take 
the form: 


H5(A;, ©) 
a {n} = >) Pet i nd 


F’ a4) e*". 
2 ie (Ai) (20) 


We will understand by the “transient response” that regimen in which there is still noticeable inherent 
motion, It is obvious that to shorten this regimen it is necessary that the absolute values of the real parts of all 
Ay be as large as possible. In the limit we shall have |Re \4| =o for all \ , i.e., the degree of system 
stability is infinite. In this case, the system's transfer function will be 


: H3(q,e) H5(q,¢) 
®, (q, 8) _ FG —_ ae , (21) 





that is, a9 = ay = ....=4]-4 = 0, and the transfer function of the system for errors, using Condition (19), will have 
the form: 
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%® )=FG 7am (22) 
Expression (22) coincides with the results in [8]. 
On the Oscillations of Sampled-Data Systems Between Sampling Times 
If y** 1, the system's forced motion, on the basis of Expression (4), is written in the form 
8 . . 
Hy (Px: ®) F, (2x) enn 
i[n, e] = >) 33 ee for O<e<y, 
k=1 G" (P%) ed (P,) sh 

(23) 
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7“ ~ Cre) ePkh, (p,) 








In the intervals Os€ < y and yse<1, the function x,[n, €] varies exponentially [6], i.e., in the general 
case the process will have an oscillatory character, But for some systems, given a definite form of forcing 
function, the system oscillations may not appear. For example, in first-order astatic sampled-data systems, 
system oscillations between sampling times are not present for unit forcing functions, since in this case H4(0,€)= 
= H7,(0, €) = const., i.e., in the steady state, the system's output function will always be constant, and not depend 
upon ¢. But for more complex forcing functions of the general form (2), Hy (Py, €) and Hn(P, » €) of Expres- 
sion (23) cannot equal one another, and will depend upon €. Consequently, in spite of the ber that Condition 
(5) for distortion-free forced motion does hold between sampling times, the output magnitude asa whole will 
have an oscillatory character, If there is a fixing element in the system, then this is equivalent to y = 1. In 
this case, $%.1(q) = $16 = $$ (q), and there will be no oscillations in the system (if the system is stable), 


Estimating the Reproduction Error of the Forcing Function 





We shall assume that the relative reproduction error of the system is determined in the following fashion: 


{ 2 [n] — kgef [nl | 
kelfl (24) 





Among the problems of designing a sampled-data system is that a limitation is placed on the system para- 
meters by the inequality 


|\e|<eo, (25) 


where €, is the given maximum error, The error in reproducing the forcing function f[n] by the forced com- 
ponent x-[n, € ] may be written in the form 


i | pln, e] | 


er[n, e] e y Rad Ve. i. (26) 


Finally, we introduce the reproduction error of the individual elementary excitations 
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2= FO” (28) 


We find from (27) and (28) 
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After some simple transformations, we find from (4), (26) and (29) 


8 F° F 
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This formula shows that to improve the quality of reproduction of the forcing function it is necessary to 
strive to reduce €( p,). 
We now turn to the investigation of:ways to reduce €(p,,). Initially we set 


e* =z. (31) 








Then, G°(p,) and H®(p,, € ) may be put in the form 
G (p,) = a, (e°* — e) (e°*# — eo)... (e*#— et) = 


i 
k 
= 4; (5, — 22,) (Zo, — 2%) ++ + (204 — 24) = aT] RE, 


- (32) 


H*(p,, ¢) = b,, (* — et) (e** — em)... (e°* — e™™) = 


= bin (Zo, — 2y,) (Zo, — 2.) « « « (2oy — Bm) = bm [] RP. (33) 
im] 


Here, rp is a Nector in the Z plane drawn from the i'th pole of the system's transfer function to the k'th 
pole of F*(q), and R, is the vector drawn from the i'th zero of the transfer function of the k'th pole of F*(q) 
(Fig. 2). Since F®(q + j2k) = F®(q) due to the periodicity [6, 7], we will consider F*(q) in the strip —4 < Lmqsr. 


The strip-« < Lmq = @ corresponds to the region lying to the right of the imaginary axis in the Z plane (Fig. 2). 
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Fig. 2. 


It is clear from Fig. 2 that rpk = RP? t RiK and poe = Re +R “*, Making these substitutions in Equations 
(32) and (33), we, obtain 


1 
G (e,) =a, |] (+R), (34) 
iw] 
HT" (9,, &) = b 11 (RP + Ri’), (35) 


i= j 


where RF, RF, and rik are the vectors shown in Fig. 2. Substituting these quantities in the expression for €(p ”) 
(27), we get 


bin I (R¥1 + Ri*) 
G* (0) i= 
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Considering that 


l 
G* (0) = a, (1 — 2%) (1—%,)...(4 —2,) = ar [] AP, (37) 


i—1 


H* (0, ©) =bm|] RF (38) 


i=—1 


and substituting (37) and (38) in (36), we find 








I m RY 
TI RP! II (rR? + Ri") Il (: a) 
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This formula allows one to estimate the error of reproduction of each elementary excitation, defined by a 
pole z py’ in the Z plane, From this it is clear that to reduce the error in reproducing the forcing function it is 
necessary that the zeroes and poles of the system's transfer function be removed from the region in which the poles 
of F*(q) are to be found, and that to the greatest extent possible the zeroes of the transfer function should be close 
to its poles, Since, in the majority of cases, the number of zeroes of the transfer function is always less than the 


number of its poles, it is expeditious to try to place the zeroes close to those poles of the transfer function which 
are closest to the region in which the poles of F*(q) lie. 
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Fig. 3 


If the region in which the poles of F*(q) lie is sufficiently small, and sufficiently remote from the region 
in toes ie poles and — of the system's transfer function lie (Fig. 3), it is possible in this case to consider 
<< 1, and Ry ak pz << 1, and Equation (39) is replaced by the simple approximate formula 


m I 
i i 
+= R"| —— | (40) 
aR Dae — 2p 
In this case, Expression (33) is also significantly simplified: 
ep[n] = y Files) SRM ePk"| , (41) 
~T7hl kai ¢* Fa(0,) 
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where 
I 


Sad —.— 5 -1., AY = (f* —1). 
Bae 2 ae 


We transform Expression (41), noting that 


8 F° 
Af [n] = f[n +41] —/f[n] = >) = 1 (Px) (e°* — 4) eh”, 
kare “Fe (p,) 
As a result, we obtain 
ss — PiAfialt 
where R, = (e °° —1), D = SRy. 
Thus, to reduce the reproduction error, it is necessary to reduce the quantity 
~ i ow 
D=R —-— >), —],. (43 
2 ine — 2 ap 


which depends only on the mutual distribution of the zeroes and poles of the system's transfer function, and on 
the quantity pg (Fig. 3). 


The Problem of Reducing the Innerent Component 








As is well-known, the inherent component is an unavoidable reproduction error during the transient response 
and, in the general case, has the form: 


I * ° 
FL (,) HQ, & 
x(n] = >; S,e™ = roy a1 Ow ei™ (44) 
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We consider initially the inherent motion due to one component of the forcing function 


fe (n] = cxe**”. (45) 


The discrete Laplace transform for this compo- 
nent will be 
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and the component of motion, x,,[n], has the expression 
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Fig. 4, 
Here, 
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a Ee = oo 


d 


fe 


fi 


fc 


fo 


in 





Thus, 


l Agn 
c,b (2,, — 2.) (2), — 2y,)---(%, —2,_) en 
rj, (n] =>) km i Y i i Ym 





49 
4m Gr, — 4a, — 2a) - > Gy — y_,) &, =A) yma) (49) 


We emphasize that here the sum extends over all poles of the system's transfer function, The expression 
obtained may be given a vector interpretation (Fig. 4). It is not difficult to see that each parentheses of the 
numerator of Expression (49) may be represented by a vector Rip drawn from a zero, zy of the transfer 
function to its i'th pole, Zi and the parentheses of the denominator of Expression (49) may be represented by 
the product of vectors RiP , drawn to the ith pole of the system's transfer function from all other poles, multi- 


plied by the vector RP , drawn from the given k'th pole of the forcing function to the i"th pole of the 
system's transfer function (Fig. 4). 


Using the foregoing notation, we can write (49) as: 
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In the case where the poles of the forcing function lie in a small region about the point (1, j0) in the Z plane 
(Fig. 3), it is possible, with sufficient accuracy, to take RP mS RP =-Rk{” (Figs. 2 and 3), In this case, Expres- 
sion (50) may be written in the form 
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Expression (51) shows that in this case to lower the inherent component of oscillation, it is necessary to 
decrease the vectors Rp, going from the poles z) of the transfer function to the point (1, j0) and from the 
zeroes z,, of the transfer function to its poles, and to increase the vectors drawn between the poles of the trans~ 
fer function and drawn from its zeroes to the point (1, j0). It is not difficult to convince oneself that these require 
ments cannot be met independently of each other and that to a certain extent they are mutually contradictory. 
However, it is possible to find a certain optimum for which the distribution of zeroes and poles of the transfer 
function permit of a minimum value for the inherent component of motion, 


SUMMARY 


1. From the technical conditions imposed on the system, the limits should be determined for the region in 
which lie the poles of the transform of the forcing function which is to be reproduced with some permissible error. 


2. The conditions for distortion-free forced motion may be formulated as follows: a) the system must be 
stable, b) the poles of the system's transfer function may not coincide with the poles of the transform of the 
forcing function, c) the spectrum of the system must be constant for all poles of the forcing function, 


If the transform of the forcing function has multiple poles, then the condition for distortion-free forced 
motion is given by Equation (12). 


3. The poles of the system's transfer function must be distant from the region in which the poles of the 
forcing function's spectrum lie. 


4, The zeroes of the transfer function are best situated close to its poles, placed as far as possible from the 
origin of the Z plane. With an increase of the distance of the transfer function's zeroes from its poles goes an 
increase in the oscillation of the inherent component, 
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5. The poles of the system's transfer function must be placed as close as possible to the origin of the 
Z plane, 
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ON IMPROVING THE TRANSIENT RESPONSE OF CORRECTING LINKS 
WITH VARIABLE PARAMETERS 


E. K. Shigin 


(Moscow) 


Some methods are presented for improving the transient responses 
in fourth-order automatic control systems containing two integrating 
links, It is shown that when a differentiator with a variable time constant 
is introduced, the transient response can effectively be improved by the 
introduction of an integrating link with a variable time constant into the 
correcting circuit. 


The behavior of autgmatic control systems can sometimes be improved by the introduction of artificial 
nonlinearities, obtained by varying parameters during the length of the transient response, This artifice allows 
one to obtain a duration of transient response which, for given regulated object and executive organ characteristics, 
cannot be attained by linear systems with constant parameters, A better approximation to an aperiodic process 
may also be obtained [1,. 5]. 


One variable parameter which may be used is the gain of a differentiator, smoothly or step-wise varying 
during the transient response [2, 3]. The introduction of one such variable-gain differentiator immediately leads 
to the desired results in systems without integrators in the 
correcting circuit, but in systems with such integrators 
the introduction of just one variable differentiator may 
not suffice to give the requisite improvement, 


To show this, we use a fourth-order automatic con- 
trol system, the block schematic of which is given in 
Fig. 1. The system can function either as a following 
system or as a stabilizing system, The numbered boxes 
in Fig. 1 are as follows: 1 is the sensing element, 2 is 
the differentiator, 3 is the integrator, 4 is the amplifier, 
5 is the executive organ and 6 is the controlled object, 
The purposes of boxes 7 and 8 will be given below. 














The link equations may be written as follows: 


Ly = Ig — Tq; (1) 

ty = (Tp + 1)%, (2) 

mp (3) 

m= K(x +my=aK(14T +q-)2 (4) 
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(7',p + 1)7% =%, (5) 

(747 5p” “+ Tsp) a = f(t) + 25. (6) 








Here, Xq is the controlling stimulus, x, is the error, x, is the differentiator output, xs is the integrator out- 
put, xg is the amplifier output, xs is the executive organ output, x¢ is the controlled quantity, T, is the dif- 
ferentiator's time constant, T, is the integrator time constant, T; is the executive organ's time constant, T, is 
the controlled object's time constant, Tg, is the time constant of the damped controlled object, K is the amplifier 
gain, f(t) is the external disturbance and p is the symbol for differentiation, 


Using standard methods [1], let us find acceptable adjustments of the differentiator and integrator with the 
condition that 0 = 20%, where o is the magnitude of the first regulating surge following a unit step controlling 
input (X) = 1 fort> 0). Here, the transient response may be either optimal, in the ordinary sense, or nonoptimal. 


The dotted curve on Fig. 2 shows the transient response of the system under consideration for Ty = 0.7 
sec, T, = 8 sec, Ts = 0.2 sec, Tg = Ts = 1 sec and K = 2, the curve being constructed by the Bashkirov method 
(with unit step excitation), With these settings, the transient time t, = 7.5 sec (we take the length of the transient 
response to be the time necessary to establish an error not exceeding 5% of the unit step excitation). 


If the differentiator time constant T, is decreased to 0.35 sec, the transient time decreases (t, < tr,) (Fig. 
2, curve A), but the magnitude of the initial overshoot o increases, 


We consider whether, in this case, it is impossible, by stepwise varying of the parameter T, during the 
transient response, to decrease the amount of overshoot without thereby increasing the duration of the transient 
response, 


The rules for introducing step-wise variation of parameter T; can be chosen according to the recommen- 
dations given in [2, 3]. However, this may more easily be done by considering the relationship, during the course 
of the transient response, of the signs of the speed and velocity errors of x; [4]. Let Ty depend in the following 
way on the relationship of the signs of the first and second derivative errors: 


| T1 for —! >O (*startup”), 
zy 


sit 1 (A) 
KasT yo for — <0 ("braking"), 
| 


where Ty is the basic differentiator time constant (Ty = 0.35 sec.) on the startup portion and Ka, is the coef- 
ficient of step-wise variation of T, along the braking portion of the system. 


The determination of the relationship between the signs of X, and X, is carried out within box 7 (Fig. 1), 
and the discriminator which stipulates the effect of A; as a function of this relationship is located in box 2, where 


If Ty is increased at time t (Fig. 2), the transient response occurs with a lower value of overshoot o (Fig. 
2, curve B), but the duration of the transient response ty, increases. 


This type of variation of the form of the transient response curve, and of the transient time, for systems 
with variable-parameter differentiators does not depend on the way in which the variations of the parameters is 
introduced, and is inherent in all systems with two or more integrators. Indeed, it is clear from Equations (1)- 
(6) that, after the new equilibrium state is established (x, = xg = 1 for t = oo), all remaining intermediate quanti- 
ties vanish: 


1 =%,=%=%=21,=0. (7) 


On the basis of Equation (3), one may write 


4 b= =] 
ah) F,\ adt = 0. (8) 
0 


For the transient response represented by Fig. 2, curve A, Expression (8) can be written in the form 
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Fig. 2. Graphs of the transient responses in an automatic control system 
with variable time constants in the differentiator and integrator. 


t, t, 
~~ i i 
a = a \nde+ 7\ nat + 7.) at + 7) "at = 0, (9) 
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where tg, tg and tg are the instants of time at which the error becomes equal to zero. 


The first and third integrals in (9) are positive and the second and fourth are negative, Since Equation (8) 
must be valid also for systems with variable T,, it is then possible to write, for the transient response represented 
by Fig. 2, curve B, 


2 co 
ff” on 5 \ aa + 7;\ z,dt = 0, (10) 
0 


, 


where tg is the moment at which the error becomes equal to zero. 
For Equation (10) to hold, the integrals must be of different sign, Thus, the transient response cannot occur 


without overshoot. 


Since the error x, in the second integral of Equation (9) is greater in absolute value than the error x, 
in the second integral of Equation (10), the value of xg in Equation (9) will approximate to its equilibrium value 
in a shorter time than will the value of xs in (10). 


It is obvious that in order to decrease the duration of the transient response at the moment ty, it is neces- 
sary to have a step-wise decrease of the integrator time constant, Tg. We stipulate that T, will vary as a function 
of the relationship of the signs of the error and its first derivative: 


Tx0 for <0, 
1 


T,= f (B) 
1 
Ky, for r > 0, 
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Fig. 3, Functional dependence of the relative duration of transient 


response and of overshoot on the variation coefficients of the differen~ 
tlator and integrator time constants, 
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Fig. 4. A is the graph of x, = 0.25 t, B is the transient response 
for a constant-parameter system, C for a variable~parameter 
system, both in the case when x, = 0.25 t, D and E are the curves 
for, respectively, constant-parameter and variable- parameter 


system transient response for a step-wise excitation, f(t) = 1, for 
t> 0. 








where Tg is the basic integrator time constant and Ka, is the coefficient of decrease of the integrator time constant 
along the interval of forced integration, 


In this case, Equation (10) gets rewritten as: 


, 


(oo) 4 Ka, 
wi \ dt - T'29 
0 





: * 
x,dt + Tx \ z,dt, (11) 


~' 


2 
where the second integral corresponds to the integration from the onset of overshoot at t, to the change in sign 
of the velocity error X; at time tg. The third integral corresponds to the integration of the negative error with 
the basic time constant Tg). Curve C of Fig. 2 is the graph of the transient response of a system with variable 


parameters T; and T, given by Conditions (A) and (B), with Ty = 0.35 sec, Ka, = 3, Tyy = 8 sec, and Ka, * 10, 
The duration t;, of the transient response here is less than t;,. 


To determine the influence of Ka, and Ka, on the character of the transient response, one may construct 
curves showing the dependence of the duration of the transient response and the magnitude of the overshoot on 
Ka, and Ka,. Figure 3 shows a certain region of these characteristics, corresponding to a transient response with 
just one essential overshoot. 


The ordinates of Fig. 3 correspond to the relative time of the transient response T = t., /t,,, where t., is 
the duration of the transient response for a system with variable parameters, and te is the duration of the transient 
response for a systern with constant parameters, The solid curves give the dependence T = F(Ky, ) for various 
values of . The dotted curves join points with equal values of overshoot. The value of this overshoot, in 
percentages, is given above each curve. The dotted curve a-b bounds the region with just one essential over- 
shoot. 


The cross-hatched region represents the zone of allowable values of Ka, and Ka, satisfying the conditions 
Oo s 20%and Ts 1, 


Real improvements in the transient responses of systems with parameters varying according to Laws (A) and 
(B) are also observed with linearly varying input quantities (x, = wt, where w is the speed of variation of the 
input quantity) and in cases of step-wise varying external excitations, when the system is working as a stabilizer. 
The corresponding curves for the transient responses are given in Fig. 4 for systems adjusted in accordance with 
point a of Fig. 3. 


SUMMARY 


1, The introduction into automatic control systems of given structure and characteristics, and with essen- 
tially invariant links, of parameters which vary during the course of the transient response can significantly 
improve the quality of the transient response, 


2. For systems with integrators, this improvement is obtained, not only by the introduction of variable 
time constants in the differentiators, but also by a mandatory introduction of variable time constants in the 
integrators, 
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STATISTICAL INVESTIGATION OF NONSTATIONARY PROCESSES IN 
LINEAR SYSTEMS BY MEANS OF INVERSE SIMULATING DEVICES 


A. V. Solodov 


(Moscow) 


A general theory is presented for the construction of 
what may be called inverse circuits, The application of 
inverse circuits to the investigation of nonstationary random 
processes by the simulation method is given, It is shown to be 
possible to couple inverse circuits with real regulators, The 
application of this methodology is illustrated by an example. 


INTRODUCTION 


The choice of parameters of a control system of given structure, minimizing the mean-square deviation 
of the total output signal from the effective signal, leads in general, when the random input signal and the 
system are nonstationary and the process is observed at time t = 0, to the problem of minimizing the expression 


‘n 


nh 
out (tn) = \ | W (tn, §)W (tn, u) Kin (u, §) dédu. (1) 
0 
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This expression defines the relationship between the dispersion 6,,,, of the output signal at a certain fixed 
moment of time t,, and the correlation function Kj, (u, € ) of the input signal, in terms of the system's impulsive 
response function W(t, € ) [1]. Due to the great difficulty of obtaining analytical solutions to this problem, in the 
majority of practical cases recourse is had to numerical solutions using the resources of modern computing tech~- 
nology. Here, simulating devices merit particular attention since they have been widely used in the design of 
practical systems, due to the simplicity of the set of problems, the graphical nature of the solutions and the pos- 
sibility of coupling with actual regulators, However, their use for statistical analysis has heretofore been limited 
to certain special cases only. If the input signal is a stationary process, of the “white” noise type, then 


th 


Kin (u, é) = Sod (t 7 ¢)and Sout (tp) one So \ w? (t, §) dé, (2) 


where S, is the noise spectral density. 


Under these conditions, the problem of determining the dispersion of the output signal reduces to the inte- 
gration of the square of the system's impulsive response function with respect to its second argument, —. Since 
the impulsive response function W(t, € ) is defined as the reaction of the system to a delta function input applied 
at time t = — , then, as obtained experimentally, W(t, € ) will contain € as a parameter, whereas it is necessary, 
in order to solve (2) by a simulating device, to have the argument t as the parameter, in order that € might be 








the simulated variable, For this, it suffices to find a differential equation with € as the independent variable 
which would be satisfied by the function W(t, — ) for a fixed value of t. It is well known that this condition is 
satisfied by the conjugate equation, By implementing its simulation for a certain fixed value of t, one can 
obtain an output signal in the form of W(t, — ) as a function of —. 


The practical implementation of this idea directly, with the aim of simulating, has still not found wide- 
spread use and is to be foundonlyina few works, although the solution of analogous problems in another domain 
of technology, the computation of ballistic correction, was carried out by D. A. Venttsel' as early as 1926 
(published in 1928 [2]. 


In 1952, the conjugate equation was used for determining, on electronic models, integrals of the form 
t 
\ | NV (t, t)| de, 
0 


connected with the solution of Bulgakov's problem concerning accumulated perturbations. Here, the technical 
details of solving the problem were worked out, as were the theoretical questions. 


Subsequently, in [3], an analogous method of simulating the conjugate equation was presented, to be used 
in evaluating integrals of type (2). This method differed from the other only in that, instead of an instrument 
for defining the amplitude of the system's reaction, it required a squaring device to be added to the simulator 
output, In [3] there was given an excellent formulation of the correct construction of models for a conjugate 
system whereby the input and output roles are interchanged for each element entering into the system design. 
For the general case of nonstationary processes, a different methcdology was presented for determining the 
dispersion at the system output, since the method described above was not directly applicable to this problem. 


Writing Expression (1) in the form 


th tn 
2. (ta) = ( W (ty, &) dé \ W (ty, u) K;,, (u, §) du 
0 0 


and introducing the notation 


I(t, §) = \ W (ty, u) Kin (u, §) du, (3) 


> 


we obtain 


tn 
afout (fn) = | W (tm )1 (tas &) ab. (4) 


We see from Expression (3) that if in the system with impulsive response function W(t, u) we introduce 
by means of some external device the function Kj (€ , t) for a fixed value of the parameter € , then the output 
signal will be the function I(t, € ). By repeating this operation for a series of values of — we can, after the 
corresponding reorganization of the results, obtain I (t, € ) as a function of €. If we then reproduce this I(t, € ) 
by functional blocks and feed it to the input of the system under investigation then, according to Expression (4), 
we can obtain the desired dispersion 0%,,; at the output. The method just described does not require the use of 
the conjugate system, and allows the investigation to be carried out with actual regulators, However, even 
leaving aside the cumbersomeness of the experiments, such a procedure does not entail the mechanization of 
the bulk of the computational process. Thus, the method of simulating the conjugate equation allows nonstatio- 
nary processes to be investigated by simulators only when the input signal is stationary, in which case the problem 
reduces to the simulation of Relationship (2). Here, the entire problem can certainly be simulated, If it is 
necessary to keep in the simulating design, as a constituent element, an actual linear regulator (the equation for 








which might be unknown a priori), then the method of simulating the conjugate equation is obviously inapplicable 
in general. For the general nonstationary case there is now known only the method of [3] which, however, is 
essentially based on the use of simulating devices only for the intermediate operations in the computation of 
Expression (1), and which does not provide complete mechanization of the analysis. 


The present work presents a method for applying simulating devices to the investigation of linear systems 
in the general nonstationary case and is based on the construction of special simulating designs which allow 
complete mechanization of the bulk of the process of computing Expression (1) and which allows coupling with 
actual linear regulators with constant parameters (for example, RC circuits), The theory of constructing such 
designs, the so-called inverse systems, is a development of the idea of the conjugate system as applied to a trans- 
formed design. 


1. Theory of Constructing Inverse Systems. 





Below we shall look for correct transformations of designs to allow us to obtain from them, as reaction to 
delta function inputs, outputs equivalent to the variation of the corresponding impulsive response function W(t, — ) 
with variable € , for a fixed value of t. Analytically, this is tantamount to the necessity of finding differential 
equations for the links, and equations for their interconnections, with € as the independent variable, said equations 
being satisfied by the function W(t, € ) which would simultaneously satisfy the original link and interconnection 
equations in terms of the variable t. Henceforth, a system with the impulsive response function W(t, & ), with 
t variable and ¢ fixed, will be called the original system and the system with impulsive response function W (t, €), 
‘with variable ¢ and fixed t , will be called the inverse system,! 


All linear systems with variable parameters may be represented as combinations of parallel and series con- 
nections, and feedback loops, of two types of links, inertial and forcing. By an “inertial link” is understood a 
linear system described by the equation 





= d‘x 
D % (t) —* = y (2), (5) 


i-o d 
where X,); is the output signal and y is the input signal, 
A forcing link is realized by differentiating the input signal in accordance with the relationship 


d? 





X(t) = >} b;(t) 


Fin 
re! a’ 


(6) 
where x;, is the input signal and x, is the output signal. 


A series connection of the stated two types of links, when x9(t) = y(t), forms a system of the general type 
described by the equation 








n ™m di 
>) a (2) Pout = >) dj (t) i (7) 
j=0 


i i 
i—0 dt dt 


Thus, the problem consists of finding rules for obtaining the design of the inverse system for the inertial and 
forcing links, connected in parallel and series, with feedback loops, of the original system, We consider an inertial 
link, The impulsive response function G(t, € ) of the original system will satisfy Equation (5) for y = 5(t —€ ) and 
the condition of physical realizability 


G(t, =| 0 for t<c&, (8) 


G(t,&) for t>&, 





1The meaning of the term “inverse system” will be clarified in the sequel, 
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where the connection between xo and y will be defined by the expression 


t 
Tour (t) = \G(t, w) y (u) du. (9) 


We apply the following technique which, further on, will also determine the rules for constructing the 
inverse system, We assume that a delta function type signal is produced at the system @utput. Based on (5) and 
(9), this gives us 


t n 
a(t—t) = (Ge, uw) 5} a (u) 8 (w — &) du. 
0 


io 


After using the well-known identity 


ae) oo 
(1) = J 1 —t)du (10) 
we finally obtain 
(A ZG (6, Han I+ ---— HIG (t, a) + 


(11) 
+ G(t, &) a, (&) = 8(t—6&). 


Equation (11), with € as the independent variable, is satisfied by the same function G(t, € ) as the original 
Equation (5). If, in Equation (11), we replace € by t,we shall have 


(—1)» 1G (t, Ban (+. — FS (t, a (O) + 
+G*(t, 8) a(t) =3E—2), o 
where 
G’ (t, &) = G(é, 2). (13) 


The role of the fixed moment of observation, tp, will be played by the quantity €. Simulation of Equation 
(12) does not give the necessary results, since the function G(t, — ), necessary for further application, will be 
determined for values t < € , while the reaction of the real system occurs only fort = €. In fact, interchanging 
& and t in Relationship (8), we have 


0 for E<c?t, 
at 
4. _ for ED>t, (14) 


which confirms what has been stated above. Since, in Equation (12), the role of t, is played by € , by introducing 
a change of variable 


which will be the current time in the model, we shall have from (12) 


Z1¢ (t, ty) Qn (tp —t)) +... + G*(t, ty) ay (tn —t) = 8(*) (15) 








and Condition (14) will take the form 


Cll 2) = | Oo ee Sh (16) 


G* (t, tn) for t> 0, 


which corresponds to the reaction of a physically realizable system. On the basis of Equations (5) and (15) it is 
easy to construct the block schematics? shown in Fig. 1,a and Fig. 1,b, Comparing these, one can easily form 
the rules for constructing the modified schematic from the block schematic of the original circuit, 


If the block schematic of the original design consists 
of integrators and adders with unit transfer coefficient, and 
of blocks of variable coefficients, then the inverse of each 
element is taken (i.e,, each element is turned in such a way 
as to reverse the roles of input and output) and, in the blocks 
of variable coefficients, instead of the argument t the 
argument t, —T is introduced, and the result will be the 
block schematic of the altered, i.e., inverse, design, 


Input 




















Equation (12) is the conjugate of the original Equation 
(5). Therefore, in those cases when the system under inves- 
tigation leads to the simulation of an equation of the type 
of (5) (or the corresponding system of equations written in 
Fig. 1. normal form), the problem reduces to the construction of 
the conjugate system, the roles for constructing which were 
formulated above [3]. 

















Consider now a forcing link. The impulsive response function of this link is defined by the expression 


Y (t, &) = bm (t) 8)” (t —&) +... + bo (t)3(t — 8) (17) 
and is a linear combination of delta functions of various orders, To make € the independent variable, we use the 
identity 


a’ 
dt” 





3 (t—%) = (—1)" 5 3(t —8), (18) 





application of which to Expression (17) gives 


Y (t, 8) = (—1)" [bm (0) 3 (¢— 8) +. - +o (1) (EB 


Interchanging € and t and introducing, as previously, the new variable € -t = t,-t = T, we will have 


¥ (5 ta) = [bm (to—*)8 (2)] +--+ bo (tn — 2) 8 (6). (19) 


Comparison of Equations (17) and (19) shows that the correct construction of the inverse of a forcing link 
remains the same as for an inertial link. 


A series connection of the links investigated gives a system of the general type described by Equation (7) 
(Fig. 2,a). The impulsive response function W(t, € ) of this system will satisfy Equation (7) for xj, = 5(t —€ ). 
On the other hand, the function W(t, € ) is determined by the function G(t, € ) from the relationship [4] 
W(t, 8) = (— 1)" SS bm) G (by 1 +. «bo GUE, ®- (20) 
2 The crosses on the connecting lines in the schematics indicate that the signals on these lines have had their 
signs reversed, 








Interchanging € and t and bearing (13) in mind, we obtain 


Ww" (t, &) =(— 1)" = (bn (1) C(t, E)] +... by (t)G"(t, 8), 
where 
W* (t, £) = WE, 2). (21) 


We now introduce the change of variable T = t, —t,finally obtaining 


W* (x, t;) = Zi [bm (tn— 2) G" (* tn)] + « «+ bo (ta —*) G*(t, ty), (22) 


and the corresponding circuit takes the form shown in Fig. 2,b. 


Comparing circuits a and b of Fig. 2, we see that the inverse system of the coupling considered is formed 
from the original system by reversing the direction of the connections, and replacing the original links by their 
inverses, 


We now consider an original system consisting of a series connection of two systems of the general type 
(Fig. 3,a). Let the impulsive response functions of these systems satisfy the equations 


nm, ™ 
» ay (t) W(t, &) = D)d;(t) 84 (t — &), 
=0 j~0 


Ps a (23) 
} s(t) Wh (t, &) = B64; (t) 8” (t — 8). 


i=-0 


Then, the impulsive response function of the series configuration, based on the integrative connection (9) (in 
which the role of the input signal will be played by the function W(t, € )), is determined by the relationship 


W (t, &) = \Ws (t, u) Wa (u, &) du. (24) 


The lower limit of the integral in (24) is affected by the condition on the impulsive response function 
W(t, t)=0 (t< 8). (25) 
We rewrite Expression (24) in the form 


5 
W (t, &) = —\W. (u, §) W(t, u) du 
i 
and carry out the rearrangement of the arguments t and €. Then, based on (21),we shall have 
t 
W(t, )= — \ Wi (t, u) Wp tw, §) du. (26) 
g 


The block schematic corresponding to Relationship (26)* is shown in Fig. 3,b. 


For parallel configurations of systems of the general type, the inverse configuration, by virtue of the 
principle of superposition, is formed by connecting in parallel the inverses of the systems forming the original 
circuit, 


We now consider a feedback loop formed by a system of the general type and an inflexible feedback con- 
nection, In this case, the correct equation is 


The minus sign in (26) disappears when the change is made to the simulation variable T = t, —t. 
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Ligy, = F (tj, + Tour); 


where L and F are differential operators. 
For xj, = 6(t—€) we shall have 


L,W (t, §) = Fé(t—&), Lb = L—F. (27) 
Equation (27) describes a certainequivalent system consisting of a series configuration of an inertial link, 


with the operator Ly, and a forcing link, with the operator F. Corresponding to the rules previously obtained, the 
inverse configuration must satisfy the relationships 


LG’ (t, §) =8(t,8), W(t, ),= Fa (t,§), L=L—F, (28) 


where the operators L} and F* are satisfied by the corresponding inverse systems, The block schematics satisfying 
Equations (27) and (28) are shown in Fig. 4,a and, in transformed form, in Fig. 4,b, In the case where a system 
of the general type is already found in the feedback path (Fig. 5,a), it is first necessary to transfer the summing 
unit to the input of the feedback path® (Fig. 5,b), after which the circuit thus obtained will appear as a combi- 
nation of the previously considered cases of a series connection with an inflexible feedback connection. 


We now carry out the construction of the inverse system using the rules studied above, obtaining the cir- 
cuit shown in Fig. 5,c, after which, by moving the branch point, we obtain the final configuration, shown in Fig. 
5,d. 


A comparison of the block schematic of the inverse circuit with that of the original circuit in all the cases 
so far studied allows us to formulate the following general rule. 


To construct the block schematic of an inverse system from that of the original system, it suffices, in the 
latter, to reverse the direction of all connections between the component systems; it is necessary to replace each 
of the component systems by its inverse, using the rules developed above, 


It is easy to show that the inverse system of the general type with constant parameters is equivalent to the 
original system, It is sufficient for this to set all the coefficients of Equations (15) and (19) constant, and then to 
apply the inversion of the transfer function to all the circuits under consideration, A consequence of this is the 
fact that links with constant parameters (for example, RC circuits or actual linear regulators) in the original 


‘Since the operator L—F (or L* —F* ) determines a link consisting of an inertial link shunted by a feedback path 
containing a forcing link [6]. 


5 Here the link Ww is the inverse of link Wp, and possesses the property \ Walt u) We" (u, &) = 8 (t —&). 
t 
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system will enter into the corresponding inverse system without any structural changes. This means that, for 
cases when there is white noise at the system input and solution of the problem leads to the use of Formula (2), 
the system equations can be a priori unknown. 


By virtue of the linearity of the systems considered, all results obtained can be extended to systems with 
several inputs and outputs, including matrix circuits, 


2. Simulating Design for Determining the Mean Square Deviation of the Full Out- 
put Signal from the Effective Signal. 








We carry out a preliminary transformation of Formula (1) which will allow it to be used for the purpose 
of simulation, The domain of integration B of the double integral in (1) is a square with side t,, given by the 
inequalities Os ust, and0=s € <t,. 


The diagonal of the square divides region B into two parts, By and By, defined by the corresponding inequali- 
ties 


0<tqusG, (29a) 
O<usti<t, (29b) 
Thus, the integral in (1) may be given in the form 
\\ou, &) dud = {\ ou, &) dude + \\ Ou, &) dudk, (30) 
B B, By 
where 
@ (u, §) = W (tp, *) W (tm 4) Kin (u, §). (31) 


It is easy to see that the surface defined by the function $(u, £) is symmetric with respect to the diagonal 
u=€ since $(u,€) = $(€, u) which follows immediately from Formula (31). Consequently, 


{\ou, &) dud = 2{\ ow, t) dude = 2\\ O(u, &) du dk. 


B B, B, 


Going over from the double integral to a repeated one and bearing inequality (29a) in mind, we obtain 


Gout (fn) = 24 W (ty §) a( W (ty, u) K,,, (u, §) du. (32) 
Setting 
H (ty, t) = \w (tn, w) Kin (u, &) du, (33) 
we shall have 
‘n 
out (fq) = 2 W (tyy §)H (try &) dé. (34) 


A comparison of Expression (33) with Formula (24) for series-connected systems of the general type with 
variable parameters allows us to make the introduction of the correlation function, Kj, (u,€) (very important 
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for what will follow), as the impulsive response function of a certain physical system with variable parameters 
since, in Formula (33), it is possible to set Kj,,(u, & ) = 0 for u< €, due to the lower limit of integration, and 
this was previously [e.g., in Formula (1)] impossible, since Ky, (u, €) #0 foru < €. 


Let 


X, (t) = Xin(t) + Min (2) 


and 


t 
Xe(t) = \ We, u) X, (u) du, 


where X;_,(t) is the effective input signal, X9(t) is the full output signal and Nj, (t) is the noise, and let it be 
required to carry out a linear transformation of the signal Xj, to the signal Xoute 1-€-» 


Xout (t) = \We (t, ) Xin (u) du, 


0 


where Wt, (t, u) determines the given linear transformation, 


Then, the quantity € = X,,,, —X» will be defined as the deviation of the full output signal from the ef- 
fective signal for the given system. 


Based on (1), the dispersion of the error, o2 , will be determined by the relationship 


th th 
02 (tn) = \ Wir (tm &) dé \ Wir (tp, u) Kin(u, §) du + 
- .* 
+ \ W (tp, e) de| (Wn, u) K, (u, §)du— 
0 0 
tn 


| Wir (tas 1) [Kin (9, &) + Kn (, :)) dn} ’ 


(35) 


where Kj,,(u, € ) is the correlation function of the effective input signal, Ky(u, € ) is the correlation function of 
the full input signal and Kgy(n,€ ) is the cross correlation function of the effective input signal and noise. 


Due to its simplicity, we carry out the construction of the simulating design for the case when there is no 
correlation between effective input signal and noise. Here, Kgn(n, €) = 0, Ky(u, €) = Kn (u, &) + Ky (u, &). 
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Bearing in mind the transformation of the domain of integration considered above, and introducing the 
functions 


t 
A, (tn, §) = (M (fa, u) Kin (u, &) du, 
g 


tn (36) 
Hy (try &) = | W (tq, ut) [K in(u, &) + Kn (w, &)] du — 2H; (ta, 8), 
0 
we can write Expression (35) in the form 
t th 
0 0 


As was shown above, Formulas (36) characterize a series connection of systems with variable parameters, 
where the correlation functions are considered as the inpulsive response functions of certain physical systems. 
Bearing this in mind, we can set up the design shown in Fig. 6,a. It is clear from (37) that the functions H, and 
Hg musthave € as the independent variable, while the circuit of Fig. 6,a must give the reactions H, and Hy as 
functions of t. To solve this problem it is necessary to construct the inverse circuit for which it suffices to use 
the rules formulated above, The circuit thus obtained is shown in Fig. 6,b, To determine the dispersion of the 
error € in accordance with Formula (37), the functions H; and Hg must be used as the input signals to systems 
with the impulsive response functions Wir (t,, & and W(tp,€). Intaoducing into these latter the simulating 
variable T we finally obtain the design given in Fig. 6,c. With an impulse at the input at time T = 0, the 
dispersion of the error will be observed at the output at time T =t,. By varying the parameters of systems W and 
W® one can so tune the system being investigated as to obtain the minimum dispersion of the error. 


























Fig. 8. 


The circuit investigated seems to be the most general. For the case of a follower system, it is sufficient 
to set Wry (t, &) = 5(t—€), 1.e., to include in the closed-loop circuit the systems Wtr and Wf, . For the case of 


a stabilizing system Wry (t,€) = 0,K;,, = 0,and the systems Wi, , Wt, and K;,, must be linked in an open-loop 
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circuit, In the case of a stationary input signal, the systems simulating the correlation functions Kip and Ky will 
be systems with constant parameters since for a stationary process K(u, € ) = K(u,—€ ), which corresponds to the 
impulsive response function of a system with constant parameters, 


The table below provides a comparison of the two methods of investigating systems with variable para- 
meters in the general, nonstationary, case. 











Method Number of times model is Number of Total number Subsidiary operations 
used in determining one variants of model 
value of dispersion investigated | usages 
Use of original k usages to determine k Introduction of the 
system to deter- points of the curve a me functions Ky, (€ , t) 
mine integral I(t, €) and I(t,, € ) from 
K(t,, €) and the external devices 
dispersion [3] 
Use of the inverse One n n None 
system 














3. Simulating the Correlation Function of a Random Process by a Linear System 





A random process may be considered as an infinite collection of random variables depending on the time 
t. Since experimental conditions permit only a finite number of observations to be made of the possible values 
of the random function, the practical use of an experimentally derived function as the correlation function only 
approximates to some degree or other the desired values. It thus becomes clear that it is unreasonable to strive 
in a simulating device for an accurate reproduction of the form of the given correlation function, but is sufficient 
merely to provide a good enough approximation, In particular, the function 


K (ty, te) = >} gi (ta) 4 (te), (38) 
i=] 


satisfying an n'th order linear differential equation with variable coefficients, may provide, depending on the 
number of terms used, an approximation to the correlation function which is as good as desired, However, in 
the majority of practical cases, it suffices to limit oneself to one or two terms of Polynomial (38). If, for 
example, 


K (ty, te) = 1 (t1) %: (ta) (39) 
must be the impulsive response function of some linear system, then 


t, 
Zour (t1) = | K (th, u) ain (u) du. (40) 
0 


Substituting (39) in (40), differentiating with respect to t, and eliminating the intermediate variable, we 
shall have 


dz 9; (41) 
om. on ; 
dt, qi (ha) Tour = K (ty, 4) Min (41) 





Substituting xj, = 5 (ty —t,) in (40) and (41),we finally obtain 


d q, (¢1) .> Pal 
a, * (ty, te) ah) K (ty, tg) = K (ty, t1) 8 (t; — ty). 





Equation (42) also determines a linear system with variable parameters which simulates the correlation 
function for values t; = ts. For stationary processes K(t;, tg) = K(t,—t,). It is easy to show that, in this case, 
the corresponding differential equation will have constant coefficients. 


Example 


In a second-order stabilizing system described by the equation 


d*z dz dz 
at (t) 5 + x, —ju + Ko®oy, = Ko® in t+ 1 in ’ 


it is required to find, by simulation, the optimal value of the coefficient Ky corresponding to the minimum 
dispersion of the output signal o? (t,) at the moment of observation t = t,. 


The input signal is a nonstationary random process with the correlation function 


Ky= S exp[— 5 (# + |. 


The given data are characterized by the following table: 





@s (t) K 8, volt” + sec, h, ger 





4 
FD 0.9 44 2 65 














According to (42), the equation of the system which simulates the correlation function takes the form 


5-Kn (t, €) + AK y (t, &) = Se" 8 (t — 8). 

The block schematics of both original systems are given in Fig. 7. In accordance with Section 3 of this 
paper the inverse circuit for simulating the stabilizing system must have the form shown in Fig. 8a. By carrying 
out the construction of this circuit in accordance with the rules studied earlier, we finally obtain the simulating 
circuit shown in Fig. 8b. By varying the value of coefficient Ky and observing the output signal one can deter- 
mine the optimal value of Ky corresponding to minimum dispersion, 


Simulation of the circuit by a type IPT-5 model gave an optimal value of K; = 1.65. A numerical solution 
of the same problem gave a value of Ky equal to 1.62, 


SUMMARY 


The use of inverse simulating devices extends the possibilities of using statistical methods for the analysis 
and synthesis of automatic control systems containing variable parameters, A great practical advantage of the 
method of inverse systems consists in the possibility it gives of complete mechanizing the bulk of the process 
of determining the dispersion of the error, and its applicability in cases where it is necessary to use, as an 
element of the simulating device, an actual linear regulator with constant parameters, the equations for which 
may not be previously known, The accuracy of the results obtained depends basically on the accuracy of the 
simulating devices used. 
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LITERATURE CITED 


[1] V. S&S. Pugachev, Theory of Random Functions and their Use in Problems of Automatic Control [in 
Russian] (Gostekhizdat, 1957), 





va 


vo 





*In 





[2] D. A. Venttsel* Computation of Variable Trajectory Elements (F, E. Dzerzhinskii Military-Technical 
Academy Press, 1928).* 


[3] J. H. Lanning and R. H. Battin, "An application of analog computers to the statistical analysis of time- 
variable networks,” Trans, IRE vol. CT-2(March 1955). 


[4] L. A. Zadeh, "The determination of the impulsive response of variable networks," Journ, Appl. Phys. 
vol, 21 (1950), 





*In Russian, 


317 








THE INFLUENCE OF FLUCTUATIONS ON THE OPERATION OF AN 
AUTOMATIC RANGE-FINDER 


I. N. Amiantov and V. I. Tikhonoy 


(Moscow) 


Errors of distance measurement arising from fluctuating 
disturbances and the problem of stability are considered for a 
system of automatic tracking using a linear approximation of 
range. 


1. Short Description of Automatic Range~-Finder Operation. 





Figure 1 provides the block schematic for a system for the automatic tracking of a target by radar range, 
a system called, for brevity, an automatic range-finder. 





























“t | Oj mcs ug 
t t 
I, 


uD ‘out. 























Fig. 1. Block schematic of the automatic range-finder. 1 is the 
time detector, 2 is the time modulator and 3 is the differential 
detector. 


The input voltage u,,, is the sum of the effective signal and the noise: uy,(t) = u(t) + €(t) (only the signal 
u is shown in Fig. 1). The effective signals are pulses reflected from the target, assumed to be stationary, and 
thus have a fixed repetition period T,. The moment the reflected pulse appears, measured from the emission 
@ the scanning pulse, is denoted by Tr. 


The time modulator develops a pair of rectangular selector pulses, the moment of appearance of which, 
during the n'th repetition Periods is shifted with respect to the scanning pulse origination T,,, by a quantity which 
depends on the dc voltage u{®) existing at the output of the differential detector during that period. The task of 
the automatic range-finder is the accurate placement of the selector pulses on the reflected pulses, The error 
signal, depending on the difference AT, = Tp ~Tp, is obtained at the output of the time detector in the form of 
a pair of pulses of different areas, 


The time detector consists of two coincidence tubes, each of which fires only upon simultaneous action 
of one of the selector pulses and the input signal. The difference of the areas of the pulses at the output of the 


318 

















uu isi 
TT Pet me 






































hit t SYS 
—_—h ———+ 
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Fig. 2. Mutual orientation of the reflected Fig. 3, Graphic determination of AT, , 4 
signal and the selector pulses, from AT. 


time detector, depending on AT,,, is measured by the differential detector and converted to » dc voltage uf? +1) 


which controls the moment at which the selector pulses appear. This dc voltage may also be used as the system's 
output signal, 


The following assumptions are made in constructing the proposed automatic range~finder model, 


1, The form of the pulses reflected from the target is approximated by a trapezoid, but the selector pulses 
were assumed rectangular, with a definite height, The feeding to one grid of a coincidence tube of a selector 
pulse sets up a completely determinate condition for the conduction of the tube: the input signal voltage, acting 
on the other grid, must exceed a certain threshold value h; to cut off the tube, this voltage must become less 


than h, The presence of such a threshold is connected with the assumption of the rectangular form of the 
selector pulses, 


2. The differential detector reacts to the difference of areas of the pulses from the time detector S, and 
Sy (Fig. 2), LiGs, 


Aup’ = ky (Sy? — SY”), (1) 


where Aut?) is the voltage increase at the differential detector output and ky is some proportionality coefficient. 
This assumption means that the differential detector does not have a firing threshold, 


3. The shift 65T,, of the selector pulses with respect to the scanning pulse in the n‘th repetition period 
is proportional to the increase in voltage at the input to the time modulator 


87, = k,Au'), (2) 


where kg is some proportionality factor. 
It is easily seen that the equation of the closed-loop range tracking system can be written as 


AT ng, = AT, —k (SS? — SM), ke = heyhey. (3) 
In the absence of noise we have, for a small displacement AT), 


sy” T sy” = 2uAT», 
and, consequently, 


AT ns = AT, — 2ku,AT, = ATa — aAT,, % = 2huy. (4) 


This equation may be interpreted graphically, as shown in Fig. 3, The quantity AT, , 4 ~AT,, is con- 
nected with AT,, via the “reflecting line,” the tangent of whose slope equals —c. For0 < a < 1 there occurs 
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an aperiodic damping of the initial deviation AT, with increasing n, For 1< a< 2, there occurs an oscillatory 
damping of the initial deviation. For a> a,, = 2 the system is unstable, since with increasing n there occurs an 
increase in the initial deviation, 


The dependence of 1s”) -s{")) on AT, represented by the reflecting line under the given assumptions, 
is actually a straight line only for small deviations AT. 


For large deviations this dependence has a nonlinear character, as a result of which the increase of the 
initial deviation in the unstable state is limited. However, the error AT,,, is so large that the system departs from 
the auto-tracking regimen, 


2. Influence of Noise. 





At the input of the system let there be, in addition to the effective signal u(t), the action of inherent 
fluctuations from the receiver output, appearing as a stationary random process, €(t). It has a correlation time 
Tc of the order of the pulse length Tp: 


4 
c= aT Tp, (5) 


where Af is the energy bandwidth of the receiver for intermediate frequencies at the level 0.5. 


The fluctuations €(t) may have an arbitrary 
4," character (with a normal, or some other, probability 
; density). However, the intensity of the fluctuations 
§ (t), characterized by the dispersion o*, is considered 
to be not too large, so that the deviation AT,, will be 
small, and therefore the possibility of a false firing of 
the coincidence tubes can be neglected. 





TE It was assumed above that the differential detector 
E measures the difference in areas sf") —s{™) of the pulses 
at the time detector output. With the presence of not 
t very large fluctuations, this difference will be a random 
SP variable, due to two causes: 1) “creeping” fluctuations at 
the vertices and faces of the pulses during the time when 
4, the coincidence tubes are firing; 2) “tremors” of the mo- 
ment of firing of the coincidence tubes due to the super- 
imposition of fluctuations on the pulse edges (Fig. 4). 
Let us first study the first cause, Let t, — Tp, t+ 
+ Tp be the firing times (the first — conducting, the 
second — cutting off) of the coincidence tubes, without 
taking into account the variations of these times under the action of fluctuations. The time interval Tp is deter- 
mined by the pulse front voltage reaching the threshold h, In the given case, the difference of areas will equal 
(Fig. 4). 
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Fig. 4. Influence of fluctuations on the functioning 
of the time detector. SP is a selector pulse. 


AS” = Sy” — S{” = 2u,AT, + O(AT,), (6) 
where 
tnttp sn hT a % 
@M(AT,) = \ $(x) dx — \ & (x) da. 
tn —AT yp +-T tn—*p 


For a small deviation, AT, << Tc, (AT) may be expanded in a Taylor series about zero, and limited 
just to the linear terms 





—_—- x, er -* 





a 
@(AT,) > (0) + oqp— Ala = 
"P 
= \ [5 (tn + 2) —§ (tn — 2)] dx + [& (tn + to) + & (tn —%)] AT a. 


We then obtain 


AS{” == QupATn + [8 (tn + %) + (tn —%)) AT n + 


' (7) 
+ (lel +2) —&(t, —2)] dz. 


We shall ignore incorrect firings (which is approximately permissible for 20 < h) and consider a tremor 
of the moment of firing of a coincidence tube due to fluctuations. Such a study (even in a first approximation) 
requires consideration of the joint probability density of the fluctuations themselves, €(t), and their derivative 
de (t)/dt [1]. However, we shall effect a simplification, Let At be the instantaneous value of the shift (tremor) 
in the moment of firing. The variation of €(t) in the interval At is proportional to At/tc. 


Since cy; = V Ae <% the function £€ (t) can be approx- 











u,& 
tang=5- “ imately considered as constant on the segment At®. We thus obtain 
P) a simple relationship for the quantity At, the meaning of which is 
} &(t,-—/] clarified by Fig. 5: 
4 TA E 
Lf > ¢t E(t,—t E(t, + 
a ton —_ At, TS { = ; ’ At, = —_—— *?) ’ (8) 


- T di lse " 
Fig. 5. Approximate determination where s = uy /Ty is the steepness of the leading pulse edge 


of the instability of the moment of The difference of areas attributable to tne shift of the firing 
firing. moment equals 
1 
ASS” = —— [8 (tn + tp) — © (tn — tp]. (9) 


The total difference will equal 
AS™ = AS + AS. 
On the basis of Equation (3) we obtain the following equation for the closed circuit: 


AT na — AT, = — 2k AT n — k [E (tn + %) —€ (tn — %) ATn + 
tp fa (10) 
+k | [E (tm —2)—E (te + 2) de + 5% [2 (tn — tp) —P (tn +B) 
Turning to a graphical interpretation of Equation (10) in accordance with Fig. 3, we can state that the 
fluctuations occasion a variation of the slope of the "reflecting line" (second term on the right side) and a shift 
of the “reflecting line” along the axis of ordinates (third and fourth terms), The variations of the slope and shift 


of the “reflecting line" are different for different n. 


In the general case, the deviation AT,, is a nonstationary random function of the discrete times t, = nT9, 
where, as will be shown below, the mean value of aT, = 0, Such a character of AT, requires that a special way 


* Here, and in what follows, a superscript bar denotes a statistical average. 


321 














be found for determining the stability, instability and errors of this automatic range-finding system. 
It will be considered that the system is stable if there is a finite limit to the sequence at: 


ox. = lim AT3. (11) 
n- @ 


The system is unstable if, on the contrary, the sequence AT? does not have a dinite limit for n-> o. 
As a quantitative measure of the error in a stable system, we take the quantity 








CAT = / lim AT? . (12) 
n- © 
8. Statistical Characteristics. 


Since the deviation y = AT,, is a step~function of the discrete times t,, * nT», Equation (10) is a stochastic 
first-order linear difference equation with variable coefficients and a forcing function [2]. We rewrite it as 





Ay + (« + BC) y = am + br, (13) 
where 
a= hus, Pra=k, b=-TE, Cin) =E(tn+%) +E(tn—*e)s 
(14) 
Mh (n) = Fee —2)— (met andes, ia (n) = & (ty — tp) — 8 (tn + tp). 


The solution of this equation is made up of the forced component y, and the inherent component yg. In a 
stable system, the inherent component is damped for n> a, Therefore we consider initially a single forced 
solution, y;(N): 


N— 
yr (N) = yr (NTo) = Dd) [am (&) + bn (A) X 
ne aoe N (15) 
x J] t—«—Be()) = % Am, 
n=k+1 m=1 
where 
m—1 
Am = |any (N — m) + by,(N —m)} [] 4—a—B0(N— J), 

~ (16) 


m=N—k, l=N—n. 


We note that A,, consists of factors with respect to moments of time separated by at least the interval T,. 
Since Tc © Tp << Tp, the individual factors are statistically independent. Therefore, in averaging over the set 
of actual values of €(t), the mean value of the product is equal to the product of the mean values of the factors, 
1.e., 





N N 
n(N)= >) An = >) lam(N—m) + by, (N—m)] x (17) 
m=) m=1 


(cont'd) 








m-—l 


x [J f—«—R(W— J). (17) 
i=l 








If — = 0, then ny (N—-m) = Ng (N—m) = 0 for any N—m and, consequently, y; (N) = 0. 
For the dispersion we obtain 


(N= % > AAj= > Ay. (18) 


Since, for i + j, 








Ai A; = [aq (N — i) + bn (N —i)) (aq, (N — 7) + by.(N —7J)) X 


i-1 ae | 


x ]] 4—«—pe(n —p)) J] 4 —«— Be (N —») =0, 








pl v=1 
then, 
ae N m—1 
yi(N) = >) fam (N — m) + by, (N —m)p [] (1—2— BN — DP. (19) 
m=] l=] 








Since € (t) is a stationary process, the quantities [a7y,(A)-+by.(A)] and [1—a—BO(A))* donot 
depend on the index A, and therefore the right member of Expression (19) appears as the sum of N terms ofa 
geometric series, with the first term c= (ay; (A) -+-by2(A)]* and ratio g = [1 — a — BC(A)}?. 








If q< 1, then a limit exists as N-> o: 


c 
i—g~ 





lim y(N) = (20) 


It is possible to show that for q < 1 the inherent component of the solution yg is damped, Therefore, in 
accordance with the previously assumed definition of the error, we have 





SAT = Vie . (21) 


The regions of stability and instability are determined, respectively, by the expressions q = 1 andq> 1, 


We now express the error in terms of the characteristics of the original process, €(t), which will be consi- 
dered, by way of an example, to be normal, with the correlation function: 


&(t,)& (4g) = oR (| t, — |). 


In the given case 


q=(1—a)®—2(1 —a) BE + BS? = (1 — a) + 


+ B® (262 + 22 (K+ %)&(A — %)) = 
= (1 — a)® + 28%o*(1 + R(2t%)] @=0), 














c= ant (A) + B88 (A) + 2abms () Ma) = 
*p *p 
=a*| | O=2) C0 Fa O— a) — EO +a) dod, + 


oO TO 











+ Bb? [& (A — cp) + EA (A + tp) — 267 (A — mp) F200 + tp) + 
‘P 
+ 20b | (PR —)—FP OF ep Ea—2)—tO +a) de = 


70 





*?*P 
= 2a*o* \ \ [R (| a — Z|) — R( |x, + 2|)[ dxydz, + 4b%o* [1 — R? (2cp)). 
0 70 
The last equality was written on the basis that, for a normal distribution with zero mean, we have the values 


B= Sot, G7 (ti) 7 (ta) = of [1 + 2R?2(|ty—ty|)], (A) E(t) = 0. 


Substituting the values for c and q in Formula (21), we obtain 











tptp 
tyr = — 11 + RW | | Rm —m)— 
0 Te Te 
20° tt (22) 
—R(\% + 2) dey dz, + SE [1 — R(2rp)| | 
0 
where 
a(2—a) (23) 





= VY BH+FRem) * 
The first term within the curly brackets in (22) characterizes the influence of creeping on the vertices and 
edges of the pulses, and the second term characterizes the shift of the moment of firing. 


As will be shown below, the stability condition q< 1 is equivalent to the condition that o* <o3, There- 
fore, in the stable regime, the quantity o,7 is real. For o = 0, Gay = 0, and for 6 = 09, 0,,, = @. For small 
© we obtain the approximate formula 


*p *p 
ohn = (<)"{t + R20) \ | LR (|x; —29|) — A (lay + 29))) day day. (24) 


In accordance with the assumptions made with respect to the smallness of the deviation AT,,, for Formulas 
(23) and (24) to be correct the dispersion o* must be sufficiently small so that the inequality og7 << Tc holds, 


Numerical estimates can be made from Formulas (23) and (24) if some concrete form of the correlation 
coefficient R(T ) is given (cf. the example given below). 
4. On the Stability of the System. 

The limits of stability are determined from the equation 


q = (4 — a)? + 2870? [1 + FR (2t9)] = 1. 





This equation has the solution 


“oy a(2—a) 
f= % = eR Otol’ 








nw oo. Oo 


~. 
} —q 





The quantity 6, characterizes the limits of stability: for 0< o9, the system is stable; for 0 > oy the 
system is unstable, 


Substituting the values of a and g from (14), we obtain the 
2 final formula for 0}; 






2 2 Uo 
= TERGw & (1 —*%) (25) 


In a three-dimensional space with the coordinates k, uy 
he and o”, Equation (25) determines a certain surface, represented 





et in Fig. 6. The intersection of this surface with the plane k = 
u, = const gives a parabola, and with the plane uy = const, a hyper- 
bola. 


Fig. 6. Region of stability with fluctua- 


tions present, The points lying inside the surface are in the region of 


stability, and those outside of it are in the region of instability. 
The line Ly on the ug, k plane, the equation for which is uy, = 

= 1/k, is the boundary of ordinary stability (o = 0, a= 2). For a fixed k, the dependence of o% on uy has a 
maximum at the point u, = 1/2k, 


Example 


Let R(t) = exp (— qe) . 


*c 


Then, carrying out the integration in Formula (22) and assuming that Ty <<TC, To << Tp, we obtain 


TE ae 
Sgt. dee] 2a. eee 
” a5 — 0° 7c Tc s 


(8 F(a) Ys 


Let the following relationships hold: T p = 2T p, Tc = Tp = 2Tp, Up =h, O = up/2, 0/09 = 0.6, We can 
convince ourselves, from estimates, that the effect of creep on the vertices and edges of the pulses is approxima- 
tely two times as large as the effect of the shift of the moment of firing of the coincidence tubes, and that the 
sum of the two effects leads to an error 


Oar ~ 0.2Tp, (27) 


If Tp = 1 microsecond, then the error in measuring distance to the target equals Op = 0 arv/2 30 meters, 
where vy is the velocity of light. We note that Equation (27) is justified by considering that we are dealing with 
small deviations, AT, << To. 


SUMMARY 


In this paper we considered the functioning of the simplest model of an automatic range~finder with the 
presence of sufficiently small fluctuations and a stationary target. Of course, there are desirable generalizations in 
the following directions: 1) rendering the model of the automatic range-finder more precise and more complex, 
2) considering the nature of target motion, 3) removing the assumptions on the smallness of the fluctuations and, 
in particular, considering incorrect firing of the coincidence tubes, 


Such investigations are necessary steps towards the making of comparative estimates of noise stability of 








automatic tracking devices by angular coordinates and distance with different forms of noise, 


Received March 29, 1957 


LITERATURE CITED 


[1] V.1. Tikhonov, "The action of small fluctuations on electronic relays," Vestnik Mosk. Univ., No. 5 
(1956). 


[2] A. O. Gel'fond, Calculus of Finite Differences (Gostekhizdat, 1952), 








DETERMINATION OF SYSTEM PARAMETERS FROM EXPERIMENTAL 
FREQUENCY CHARACTERISTICS® 


A. A. Kardashov and L. V. Karniushin 
(Lvov) 


A method is given for determining the parameters 
of linear links and automatic control systems from their 
approximate experimental amplitude-phase characteristics 
(a.p.c.). By interpolation, approximate values are found 
for the values of the coefficients in the analytic expression 
for the a.p.c., and then corrections to the coefficients thus 
found are computed by the method of least squares. 


The method to be given is applicable to transfer 
functions with arbitrary polynomials in numerator and 
denominator, and gives sufficiently accurate results 
when calculations are done by slide rule. 


Examples clarify the computational steps. 


In the literature on the theory of automatic control [1-3], there are investigated only the simplest cases 
of approximate determination of the parameters of automatic control system (a.c.s.) links from empirically~ 
determined dynamic characteristics, In recent years, a series of new works, detailing further development of 
methods for solving this problem, have been published ([4-6] and others), 


The present paper attempts to give a universal and sufficiently accurate method of finding the values of 
the coefficients of the transfer function, or of individual parameters, in linear models of real links and a.c.s., 
using approximations to their empirically determined a.p.c. 


STATEMENT OF THE PROBLEM 


To determine the desired parameters of links or of link circuits in automatic control systems, one may 
start either from given amplitude-phase characteristics or from experimentally-determined amplitude-phase 


characteristics, obtained under conditions close to the actual ones. and smoothed by well-known statistical 
means. 


An approximate structure of the amplitude-phase characteristic may, in the majority of cases, be deter- 
mined beforehand on the basis of known physical laws to which the functioning of the links (link circuits) of 
the automatic control system** are subject. In case these laws are not sufficiently well known, the necessity 
arises of determining an acceptable structure for the analytic expression of the amplitude-phase characteristic. 
This necessity may also arisen when the link equations have been simplified as a preliminary step. 


* The basic statement of this work was made on April 11, 1955, at the Fourteenth Conference of the Lvov Poly- 
technic Institute on the results of research work carried out in the year 1954, 

* * The experimentally-determined characteristic in the time domain can also be used, it being a simple matter 
to construct the amplitude-phase characteristic from this, 








If the physical relationships and initial conditions allow only the structure of the numerator of the link's 
transfer function to be determined, then the degree n of the denominator may be found from the experimental 
amplitude-phase characteristic, using the argument principle: 


n=1+m’—m", (1) 


where ! is the number of quadrants traversed by the amplitude-phase characteristic vector when the frequency 
w varies from 0 to oo, m’ is the number of roots of the numerator with negative real parts and m" is the number 
of roots of the numerator with positive real parts. A root of zero in the numerator leads to a counterclockwise 
rotation of the amplitude-phase characteristic vector by the angle » =5 m***, where m™ is the multiplicity of 
the zero root. 


In the general case, the structure of the initial linear differential equation will have the form 


(anp™ + Qn—.p™! +++-+ ap + 1) Your = 


(2) 
= (bmp™ + bm prt + 93+ bp + k) Yin» 


where m< n, Ujp and u,,,, are the input and output coordinates, a,, b, (s = 1,2,...,n;r = 1,2,...,m) andk 
are constant coefficients, and p = d/dt is the operator for differentiation with respect to time. 


To Equation (2) there corresponds the amplitude- phase characteristic 


K(jo) —- ¥n(®) + 72, () 


W (ju) = D(jo) ~ “uge) + jvy(@) ~ U (a) + Vo), ©) 





where 


Un (@) = k — dyes? + byeast— bgw® +--+, Vn (w) = bo — dyw* + byw'— byu?+---, 
Ug (@) = 1 — ages? + ayeos* — agn® +--+, vg (@) = ao — aye + age’ — a,0"- +, 


The problem reduces to finding those numerical values of the coefficients as, b, and k and, consequently, 
the desired link parameters, for which Function (3) will coincide as closely as possible with the experimental 
amplitude-phase characteristics. 


Connected with the linearization process, and also with the realization that, in practice, it is not possible 
to take into account all the physical phenomena bearing on the transient responses of automatic control system 
links, it must be borne in mind that the function, determined from analytic Expression (3), with its given structure, 
may not coincide with the experimental amplitude-phase characteristic even with selected values of the coef- 
ficients as, by and k, It follows from this that the problem of determining these coefficients may be solved by 
simple interpolation only in individual cases, and that, therefore, in the solution of the problem under consider- 
ation, use must be made of a method giving the best, in some sense or other, approximation to the function. 


For a link with a static characteristic, the coefficient k is the gain coefficient, which is most simply 
determined from the initial portion of the amplitude-phase characteristic: 


k = W (ja) | o-o = U (0) (4) 
or from the characteristic in the time domain 


u_.., (co) 
wire Uj np (90} ? (5) 


For a link with an astatic characteristic, i.e., when D(p) = ap" +8,-3P"" 1. ...+a,p, the coefficient 
k may be found by the same method, to be given below, as the other coefficients. 








If the value of k is found by a preliminary step, then the solution of the problem entails setting up and 
solving n +m linear algebraic equations in the n + m desired coefficients a, and by (n + m + 1 equations for 
an astatic system). These equations can be set up using known principles of function approximation [7]. 


DESCRIPTION OF THE METHOD 


One may present as follows the nub of the method for determining the numerical values for the desired 
coefficients ag and b;: initially, rough values of these coefficients are determined using the method, detailed 
below, of simple interpolation from the experimental (given) amplitude-phase characteristic, expressed in linear- 
fractional form (3), and thereafter, if necessary, the approximation to the amplitude-phase characteristic is 
brought into better correspondence with the experimental curve at intermediate points by improving the coef- 
ficients already found by means of least-squared computations, 


Method of Interpolation 





In the interpolation method given here, the interpolation points are chosen at arbitrary points, not at the 
intersections of the amplitude-phase characteristic with the axes of the phase plane, as is the case in the method 
given in[5]}. Since the interpolation is performed only as a preliminary step in the determination of the desired 
coefficients, the choice of interpolation points does not exert much influence on the accuracy of the final result. 
Still, it is advisable that the chosen points lie within the relevant frequency band, 


To find numerical values of the n+ m coefficients a; and b, it suffices to choose N points on the experi- 
mental amplitude-phase characteristic, where N = ; (n + m) forn+m even, and N = H (n + m + 1) for odd n+ m, 
Starting from the condition 
W (jo) — We (jw) =0, (8) 


where W,( jw) = Ue(w) + jV_(w) is the experimental amplitude-phase characteristic, and taking (3) into account, 
we obtain the following system of equations for the N interpolation points: 


Un(@;) + [2% (@;) 
ug (@;) + 1% (@;) 





—U, (@;) — JY («;) =0 (i = 1, 2, acs’ N). (7) 


After multiplication by the common denominator, and separation of real and imaginary parts, System (7) 
takes the form: 


Uy (01) Ug (i) — vg (i) Ve (@i) — Uy (i) = 0, 


- (8 
u g(r) Ve (601) + 04 01) Ug (ci) — Dg (ea) = 0 4,2, sry 


Equations (8) are linear, and there are 2N = n+ m of them, i.e., there are as many equations as there are 
desired coefficients.* The degree of the equations of System (8) with respect to a, and b, may easily be reduced 
by one. For this, we first eliminate vj(w,) from the first n equations, and then eliminate ug(u), obtaining 
the following two systems of equations: 


ug (aj) (U; (oj) ot. - («;)] — Vn («;) A («;) — un (@;) U, (o;) =Q (i=1. 2, vey n), (9) 
v, (ca) (U2 (ci) ++ VE (c))] 04 (01) U, (5) + ty (04) Ve (Cs) = O12 00m (0) 


the solutions of which give the dependence of the sought quantities, ay, ag, ag,..., and also ag, ag, ag... ON the 
quantities by, by, bs,...,. We then substitute these relationships in the remaining m equations of System (8), 





* When n + m is an odd number, it is necessary to choose just one of the two equations of System (8) for some one 
of the chosen points. 








obtaining a system of equations with respect to the unknowns by, bg, bs,..., the solutions of which will be the 


numbers bi. by, bg... . By means of these, one easily determines the numerical values of a}, ag, a$,... eae 


After the numerical values of the desired coefficients are found, it is possible to construct the computed 
amplitude-phase characteristic for comparison with the experimental. 


To the extent that the method of simple interpolation uses a limited number of interpolation points, coin- 
cidence of the experimental and computed amplitude~-phase characteristics may be expected only at the chosen 
points, and significant discrepancies may be expected at the intermediate points. 


Therefore, in those cases when the coefficients found by the method of simple interpolation do not provide 
a sufficiently accurate mirroring of the dynamic properties of the links, it is necessary to change over from the 
interpolation method to a method for approximating the experimental amplitude-phase characteristic by an 
expression for this characteristic which corresponds to the given differential equation. 


Computed Corrections 





To obtain improved coincidence of the experimental, W,(jw), and computed, W(jw), amplitude-phase 
characteristics, not only at the interpolation points but also at intermediate points, it is necessary to add cor- 
rections, Aa;, Aag,..., Aa,, Aby, Abg,...Ab,,, to the numerical values of the coefficients as found by the 
method of simple interpolation, These corrections are most expeditiously found by the method of least squares 
{7, 8], i.e., from the condition: 


! | K (jo) 
py Dae, ~ “ 





(i>), (11) 





where J is comprised, not only of the previously used points of the experimental amplitude-phase characteristic, 
but also of some intermediate points. 


Condition (11) can be put in the following form: 


t 
>; (0 @i) — U, (a) }* + (V (ai) — Ve (@;)]?} = min. (12) 


i=l 


If we assume that the desired corrections are small, we can linearize the functions U(w,;) and V(44), 
expanding them in Taylor series and restricting ourselves to terms of the first power in the quantities Aa;, Aag,..., 
Map, &b;, Aby,...,Ab,,. Then, instead of Expression (12), we obtain 


l m - OU (we i) aU (a) 


» {> ee. Ab, + > 5a, ) Aa, — AU («)) 4 


i—1 r=] 


+[3 (Te) Ab, +D (Tied. ) Aa, — AV (a) |} = mis, 





(13) 





where 
AU («xi) = Ug (wx) — Ug (wr), AV (0) = Ve (ax) — Vo (on). 
From known rules for minimizing functions, Condition (13) transforms to the system of equations: 


{LS ar”), am + 33 (Man) 80 — a0 con] Ma”), + 


r=1 











(g = 1,2, ..., m), (14) 


+[3 Cae *) Ab, + 3 3 (ae) 40 .— AV (o (5 =) | =0, (cont'd) 


rT=1 ~a4 








ie: 











SSCP FC an 
eae o (h==1,2,...,2 4) 


+13 Coe), ah + SY) am, ao] SSP) wo 


OV (;) dV (@,) 








Using the fact that, for example, 


OW (ja;) +t: aU («,;) . V (a) 








0b, Ob; ‘a (15) 
we may write System (14) in the form 
: / OW (ja,) n dW (ja,) ’ dW (ja,) 
Re [53 (Gg) Abe 3S (Mae), ee — AW Gio] (Bem) } =O 
1 r= s=1 


(g =1,2,..., mm), 














(16) 
aW (ja,) OW a - ) aw wos 
3 mo (E) ab + 33 CM) na, a0 as] MOP) mo 
= Tt =] 
(hk = 1,2,..., M). 
Making use of (3), we find 
OW (ja,) (ja, )" 0 (4) if cen VP 
( i i ) on De tia) = B, (ca,) (a) (ja,)’, 
(17) 


dW (je Ko a ieee 
aE) ee en (jai)* = Ap (wi) e™™* (je)’. 


We substitute the values of the derivatives from (17) in the system of Equations (16), and also represent 
4W(jw;) in the form 


AW (jen) = Co (cx) ef, 


obtaining 


l 


3 Re {| By (a) ef" > (jeo)" Ab, + Ag (co) eo? 3 (ja,)* Aa, — 


i=_] 


i Re (ca) eibe (@;) B, (ca) ore (— jon)*} =—0 (g = 4, 2, ..., m), (18) 


y Re (|B. (c) ei (@;) 3 (ja) Ab, Y Ag («) ei” (@;) (jen,)* Aa,— 
i=l r=1 s=1 


cnties Cy («;) eihe i] Ay («;) Pr eet (— jen)" =0 (h as ‘, 2, ijt n)e 


We multiply out the parentheses and brackets and, for simplicity of writing, we omit the notation for the 


ranges of summation and for the indices of the coefficients, finally obtaining the system of n + m normal equations 
in the desired corrections: 








>) Brw*Ad, — >) BrwtAb, +>) Bw Ab,—- - «+. > AB cos (2 — 9) w*Aa, — 
— >) AB sin (a — 9) w*Aa,— >) AB cos (a — 9) w'Aa, +>) AB sin (2—o) w®Aa, + 
-+- >) AB cos (a — 9) w®Aa, —---- = > CB sin (B — 9) a, 
>) B*wtAb, — >) B*w*Ab, + >) B*w8Ab, —-- - + >) AB sin (« — ) w®Aa, + 
+>) AB cos (« — 9) w*Aa, —>) AB sin («— 9) w®Aa, — >) ABcos (a — p) w®Aa,+ 
+ >} AB sin (a — 9) w"Aa, + --- = — 3} CBcos (8 — 9) ®, 

— > BrwAb, + >} B*w*Ad; — >) B® Ab, + ---— 5) AB cos (a — 9) Aa, + 
+) AB sin (« — ¢) w®Aa, + >) AB cos («—¢) o*Aa, — >) AB sin(« —¢) w’Aa,— 
— >) AB cos (a — ¢) w®Aa, +---= — >) CBsin (B — 9) 0°, 

— }) B*a*Ab, + >) B*a*Ab, — >) BYo Ad, + ---— >) AB sin (2 — 9) wAa, — 
—)}} AB cos (a — 9) w*Aa, + >) AB sin (a— ¢) «Aa, + >) ABcos (a — ¢) w®*Aa,— 
— >) ABsin (« — 9) w*Aa, —---= >) CB cos (B — ¢) wt, 
>) Ba Ab, — >) B%w*Ad, + >) B*w!Ab, —---+ >) AB cos (a — 9) w*Aa, — 
— >) AB sin (a — 9) w’Aa, — >) AB cos (a —¢) w*Aa;+ >} ABsin («— ¢) w*Aa,+ 


+ > AB cos (a — ¢) oAa, —--- = >) CBsin (8 — 9) a, 
oe ee - tf S & es es € © 8 e@e@esee8es 82 6 ¢ 80 @ 2&2 6 6 Se *eeeeee*e (19) 


> AB cos (a — ¢) w*Ab, + >) AB sin (« — ¢) w®Ab, — >) ABcos (a — 9) w'Ab, — 
— >) AB sin (a — 9) w®Ab, + >) AB cos (« — ) w®Ab, +--+ >) A®w*Aa, — 
— >} A*w*Aa, + >) A2w*Aa, —---= >} CAsin(B — a) w, 

— >) AB sin (2 — 9) w®Ab, + >) ABcos (a— 9) wtAdb, + >) ABsin («x —¢) w®Ab,— 
— >} AB cos (a — ¢) w*Ab, — >) AB sin (x — 9) w7Ab, +--+ + >) A®w*tAa, — 
— »} A*w*Aa, + >) A*a*Aa,—--- = — >) CAcos (8 — a) w?, 

—) AB cos (a — 9) w*Ab,— >) AB sin (a — 9) w®Ab,+ >| AB cos(a— ¢) w®Abs+ 
+>) AB sin (a — 9) w7Ab, — >) AB cos (« — 9) w®Ab, —- --—)) A*w*Aa, + 
+>) A%w*Aa, —>) A*a*Aa, + --- = —)} CAsin (8 —a) 0, 
> AB sin (a — 9) w®Ab, — >) AB cos (« — ¢) w®Ab, —)>) AB sin (a — 9) w’Abs + 
+>) AB cos(a— )w*Ab, + >) AB sin (a — 9) w® Ab, —- --—)) A’w*Aa, + 
+ >) A%a®Aa, — >) A%a! Aa, + - -- = >) CAcos (8 — a) w', 
> AB cos (a — ¢) w*Ab, + >) AB sin (a — ¢) w7Ab, — >} AB cos (« — 9) w®Abs — 
—)> AB sin (« — ¢) «Ab, + >) AB cos (a — ¢) @ Ab, + +--+ >) A?w*Aa, — 
— >) A%a* Aa, + >) A2w”Aa, —--- = > CAsin (8 — «) @, 


Insofar as System (19) was set up on the basis of certain assumptions, since the functions U(w;) and V(w;) 
were linearized in order to obtain (19) from (12), the process of finding the corrections (if required by accuracy 
conditions on the’ coefficients of the differential equation) may be repeated, taking into account the corrections 
already found, 








Computation of the coefficients of the equations in System (19) by means of a slide rule [9] is quite simply 
accomplished, 


To render clearer the manner in which the method presented here is applied in practice, and also to evaluate 
the results thereby obtained, examples are given in the Appendix of the evaluation of the coefficients of the dif- 
ferential equations for certain links in systems, the transfer functions for which are so constituted that the inter- 
polation method expounded in [5] is inapplicable. 


SUMMARY 


The method presented here for determining the parameters of linear links and automatic control systems 
from the experimental amplitude-phase characteristic appears more universal, and more accurate, in comparison 
with other methods of solving this problem. 


As is clear from the examples given in the Appendix, in those cases when the desired parameters are deter- 
mined with sufficient accuracy using simple interpolation from the experimental amplitude-phase characteristic, 
the interpolation method expounded in this article is preferable to the method given in [5]: the method proposed 
heze allows solutions to be found for a broader class of approximating transfer functions and, moreover, does so 
without significantly complicating the necessary computation, 


A significant increase in the amount of computational work is required only in those cases in which the 
accuracy conditions on the approximate function require’ that corrections be computed for the desired coefficients. 
However, even this work, in the majority of cases, does not introduce significant toil if the computing scheme given 
in this article is employed. 


In the simplest cases, the method given here may be inferior to other methods by virtue of the amount of 
computational work required. 


The method may prove to be useful for verifying the admissibility of simplifications of the equations for 
objects and elements of automatic control systems, and for determining the equations and parameters of auto- 
matic control systems from their given dynamic characteristics, 


APPENDIX 





Example 1 





We find the coefficients of the differential equation of a link consisting of electronic and electromagnetic 
amplifiers (EA~EMA), using an oscillographic determination of the transient response (curve a on Fig. 1). 
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Fig. 1. Transient response of the EA-EMA link: a)is the experimental curve, b) 
is the computed curve. 


The circuit diagram is given in Fig. 2.° 





* The tests were carried out on an actual paper-board-making machine of the Zhidachevski paper-board combine. 





The differential equation for the link under consi- 















































Z deration has the form 
i 
ema T4T, Tat+T, 
(3 ops. 4. — p+i\i= 
+ Bhea*ema + Blea EMA 
BUA ea KEA 
EA =(TqP + 1) iF Skeateua Me 
CR $ 

+ Yin where Tg is the time constant of the longitudinal circuit 

fue aX of EMA, T,, is the time constant of the transverse circuit 





of EMA, 6 is the feedback rigidity coefficient, kp, is the 
transfer coefficient of the electronic amplifier (equal, in 
Fig. 2. EMA is the electromagnetic amplifier, EA the steady state, to i/ug), kema {8 the transfer coefficient 


is the electronic amplifier, CR are contacts of a of the electromagnetic amplifier (equal, in the steady 
relay, Py and Py are oscillograph pickups, ma is a state, to Up,,4/i). 
milliammeter. 


The approximating differential equation has the 
form: 


(agp + ayp + 1) i= (byp +k) u., 


Such a structure for the equation accords completely with the character of the amplitude-phase character- 
istic, constructed from the experimental curve of the transient response (Fig. 3, curve a). 














Fig. 3. Amplitude-phase characteristic of the EA-EMA link: 
curve a is constructed from the experimental curve of the transient 


response, b is: found by means of interpolation, curve c incorporates ' 
corrections found by the method of least squares. 


The coefficient k is determined directly from the transient response curve (Fig. 1): 


k= i. 0.31 ma/volt. 
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To determine the remaining coefficients of the differential equation, we choose two points on the amplitude- 
phase characteristic, w = 50 and w = 300, and set up, in accord with (9) and (10), the system of equations: 


(14 — 25x 10%a-) 5.30 = 0.0535 x 50 b, + 2.30%0.31, 


5.30 x 50 a, = 2.3050 b, — 0.0535 x 0.31, 


(1 — 9x 10* a.) 0.087 +- 0.616 x 300 a, — 0.31 = 0. 


Solving the given system, we get: ag = 0.000339, a} = 0.0161, bj = 0.036, 


The amplitude-phase characteristic corresponding to these coefficients is shown on Fig. 3 (curve b). 


In order to obtain a better correspondence between the experimental and computed amplitude-phase 


characteristics at points with equal w, it is necessary to determine corrections 


Aa,, Aa; and Ab,. For this, we 


choose three points on the experimental amplitude- phase characteristic: w = 50, w = 150 and w = 300. The 
computations which follow, using (19), are given in Tables 1 and 2. 


The system of normal equations is set up from the data of Tables 1 and 2: 


4280 Ab; —8400 Aa; —53,600 Aa, = 13, 8400 Ab, + 18,800 Aa, = 0.84, 
—53,600 Ab, + 5330 x 10 Aa, =—1570. 


Solving this system, we find that Aa, = —0.000004, Aa, = 0.0113 and Ab, = 0.0254. 
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Fig. 4. G is the generator, E is the exciter, EMA 

is the electromagnetic amplifier, ST is the stabilizing 
transformer, CR are relay contacts, Py and P, are 
oscillograph pickups, Vy and Vg are voltmeters, 


Carrying out a second approximating procedure, 
taking into account the corrections already found, we 
obtain Aag = 0.672 x 1075, Aa; = 0.294 x 1078 and 
Ab; = 0.206 x 10°, Hence, the desired coefficients 
of the differential equation are: ag = 0.000342, 
ay = 0.0277 and by = 0.0635. The differential equation 
itself takes the form: 


(0.000342 p* + 0.0277 p + 1) i = (0.0635 + 
+ 0,31) Ups 
The solution of this for a step-function excitation, 
U,= 11.6 volts, will be 


i=3.6+56e “ cos (35.9 t + 1.64). 


The corresponding amplitude-phase characteristic 
and transient response curve are shown on Fig. 3 (curve 
c) and Fig. 1 (curve b). 


Since the corrections found during the second 
approximation procedure are significantly less (in 
absolute magnitude) than the coefficients of the dif- 


ferential equation found by the interpolation method, and also less than the corrections found in the first approxi- 
mation procedure, it might not have been necessary to find them. 


From the coefficients of the differential equation thus found, it is easy to determine the individual para- 


meters of the EA-EMA link. 


To do this, we set up the system of equations; 








Fat 0.000342 
1+Bkgakema 
T kp 
q EA 
— 0.0635, 
1+ Bkea*emMa 








T,4+T 
: Sl oe 0.0277, 
1+ Bkpakema 
k 
EA = 0.31, 
1+ Bhea*ema 








the solution of which is* Tg = 0.0131 sec., Ty = 0,205 sec., kp, = 2.44 ma /volt, kema = 15.5 volts/ma, 


Example 2 





We consider a system for automatic control of a voltage generator, the system circuit being given in Fig, 4. 


We will determine the coefficients of the differential equation for this system by means of the transient 
response curve (Fig. 5, curve a) taken off by oscillograph. 

















rT, u t, sec 


Fig. 5. Transient response in system for automatic control of a 


voltage generator: a is the experimental curve, b is the compu- 
ted curve, 


The generator G and the exciter E may be con- 
w-5 web? ‘ m sidered as aperiodic first-order links, the electromagnetic 
L amplifier as an aperiodic second-order link, and the 
stabilizing transformer ST may be characterized by a 
transfer function of the form 








Tip 
Tip +1° 


We consider the differential equation of the system 


(asp* + agp! + agp + agp* + ap +1)ug = 





tee : = (bip +k) ue, (A) 
Fig. 6. Amplitude-phase characteristic of the The term bp on the right side of the equation stems 
system for automatic control of a voltage from the-presence of the stabilizing transformer in the 
generator: a is the transient response curve system, 
constructed from empirical data, b is the curve 
= litude~- character- 
found by interpolation. From the character of the amplitude~phase 


istic, constructed from the empirical curve of the transient 
response (Fig. 6, curve a), it follows from (1) that, for the 
given right member, the left side of Equation (A) must be, not fifth-order, but sixth-order. Since the coordinate 
values of the amplitude-phase characteristic are relatively small in absolute value in the fifth and sixth quadrants, 
the system may be described with sufficient accuracy by differential Equation (A). 


The gain coefficient k is determined from the transient response curve: 


k — “gio) == 5.35. 
Ue 


* The feedback rigidity coefficient, g = 0.182, was previously found experimentally, 








To determine the remaining coefficients of the differential equation, we choose three points, w = 2, 
w = 4 and w = 6, on the experimental amplitude-phase characteristic, and set up the following system of 
equations, making use of (9) and (10): 


(1 — 4a, + 16a,) 39,2 = — 4.55 2b, + 4,30X5.35, 
(1 — 16a, + 256a,) 154 = — 8.55 4b; — 9,005.35, 

(1 —360a_ + 1300a,) 24.0 = 4.441 6b, — 2.41X5.35, 
(2a, — 8a3 + 32a5) 39.2 = 4.30 2b, + 4.55 5.35, 


(4a, — 64ay + 1020a,) 154 = — 9,00 4b, + 8.55x5.35, 
(6a, — 216as + 7770a;) 24.0 = — 2.41 6b, — 4.41X5.35. 


The solution of this system is: ag = 0.000336, ag = 0.00254, ag = 0.0293, ag = 0.126, ay = 0.445 and b, = 
= 0,210, 


The desired differential equation can then be written as 
(0.000336 p* +- 0.00254 p* + 0.0293 p+ 0.126 p* + 0.445 p + 1) ug = (0,210 p + 5.35) ue. (B) 


The amplitude~phase characteristic corresponding to this equation is given in Fig. 6 (curve b), 


Since there is satisfactory correspondence between the experimental and computed amplitude- phase 
characteristics, it is not necessary to carry out the computation of the corrections, 


Integrating Equation (B) with a step-function excitation, u, = 24.3 volts, we get the equation for the 
transient response 


ug = 5,35 — 3.31 e995! 4 5,70 e—0-898 ! cos (4.07 t — 4.31) + 1.16 e—**™ cos (7.07 — 1.37), 


which is constructed on Fig. 5 (curve b). 
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OPTIMAL FREQUENCY DEVIATION IN ONE-CHANNEL 
TELEMETERING SYSTEMS 


Iu. 1. Chugin 


(Moscow) 


A method for calculating the optimum frequency deviation in one- 
channel telemetering systems with fluctuating interference is presented, 
based on analysis of the noise energy spectrum, 


1. Introduction 





One of the principal indications of the quality of a telemetering system is its accuracy of measurement. 
Errors in the transmission of information may arise either due to instability of coefficients in the transmitting 
element of the system or due to the influence of noise. It is well known that the use of frequency modulation 
allows a reduction in errors at the cost of increasing receiver band-width. However it is impossible by such 
means to improve without limiting the system's stability under noise or interference. At a certain value of band- 
width (or frequency deviation), the stability under interference deteriorates, Until now, no theoretic means existed 
for finding the optimal value of the frequency deviation as a function of the parameters of the signal and of the 
channel noise, thus rendering difficult the computation and construction of such systems, 


The present paper shows that, for telemetering systems with idealized characteristics of the receiver, it is 
possible to obtain an analytic expression for the magnitude of the optimal deviation, The methods of the theory 
of random processes allow this problem to be solved. 


In the analysis, the criterion used to measure the interference stability of a telemetering system is the mean- 
square error engendered, determined by the ratio of the effective noise voltage at the receiver output, in the band 
from 0 to F, to the maximum voltage of the output signal, 


2. Frequency Spectrum at the Frequency Receiver's Output 








We consider the receiving devices of a frequency receiver (Fig. 1), consisting of an input filter , which 
passes the band fy + Af, limiter 1, frequency discriminator 2 and ideal output filter @ for the receiving band 
from 0 to Fr, i. e€., the receiver consists of two linear systems (filters @; and 3) and two nonlinear ones (limiter 
and discriminator). Ordinarily, for telemetering systems, Fr <<+ Af. 


It will be assumed that the nonlinear systems possess no inertia, that the amplitude limiter is ideal, and the 
frequency discriminator is such that its output voltage is directly proportional to the deviation of the frequency 
from the central value, wy: 


Bout = (@ — Wp). (1) 





Input Output 
oe 4 aa / ka 2 q, — 
4 t4f OF ¢ 





























Fig. 1. Block schematic of the frequency receiver. 
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In order to determine the error at the receiver output, it is necessary to find the noise power in the band 
0 toF¢ of the output filter, for which, in turn, it is necessary to know the energy spectrum of the noise at the 
output of the frequency discriminator. The energy spectrum can be obtained from the Fourier transform of the 
auto-correlation function of the discriminator output voltage. 


In [1, 2] were obtained general expressions for the correlation function at the putput of the discriminator 
when the elements of the receiving devices had the characteristics given above, and when the signal frequency 
Ws, = Wy. 


If the frequency characteristic C(w) of the filter input is Gaussian, 


C (w) = C, exp |-- “ae |: (2) 


where 4w = 2nAf and 2Af is the receiver band-width, the expression for the correlation function for the noise 
R(T) may be written in the following way [2]: 








B (x) = B, (t) + By (e). (3) 
—. z 
B,(t) =a [ars + e #1 —2uc?-4) £, (+), 
B,(*) =« {ect (— as + tte 5 iH) + 
+e F(a) [a(t — ants ° 


z= (U;,/U,¥ is the ratio of signal power to mean noise power at the output of filter 4, (at the limiter's input), 
p is the correlation coefficient of the noise at the output of filter #4, E,(x) is the integral exponential function 
¢abulated in[3), p = p(T) =e", a =n (2Aff, 


The energy spectrum of the noise, W(f), is determined from the well-known relationship 


W if) =4 \ B(x) cos 2nfrdr. (6) 





Figure 2 gives the curves for the noise energy spectrum at the discriminator output for Gaussian frequency 
characteristic of the filters and with w, = w) [2, 4]. The spectrum wascomputed from Formula (6) using numerical 
integration, With a deviation of the frequency (ws =~ wp), the intensity of the noise energy spectrum, in the region 
f/Af < 1 is lowered [5], and the interference stability of the system increases, Below we shall investigate the 
worst case, when W, = we. For telemetering systems where Fs << Af one can obtain analytic expressions for the 
noise energy spectrum and power at the receiver output. Figure 2,b shows an enlarged view of the curve of the 
spectrum W(f) for z= 5, It is clear from the figure that for z = 5 (and z > 5), thespectrum, determined by 
Formula (6), is conveniently represented in the form 


W(f)=W(0)+W(/)*, (7) 


Ww (0) =4\ B(x) de. 


We write (6) in the form 




















W (1) =Wi() + Wl), 








where 
W, (f) =W, (0) + Wi (f)" = 4\ By (2) cos nied, oy 
0 
Wa(f) = Wa (0) + Walp)" = 4\ Bs (2) cos 2nfede. 
0 (10) 
Then, 
W (0) = W,(0)+W2(0), Wi) =W,(/) + Walp: (11) 
Determining the spectrum W,(f) from Formula (9), we obtain for z = 5 (Cf. Appendix I) 
3 
_ dnt Aj /f \2 (n — 1)! in/f2 
Wi(f) == (x) > a exp[— = F(a) |- (12) 


n=—1 


It is obvious from (12) that W,(0) = 0; hence W(0) = W,(0). 


Comparing the two spectra from (11), Wy(f)* and W,(f)* , one may show that W,(f)* << W,(f)*. To find 
these spectra, we multiply the functions By(T) and B,(T ) by (cos 2nfT —1), and integrate the resulting expressions, 


wit) 
fs or 2At 





For function Bg(T ), only numerical integration is possible, 
The computations confirm that the spectrum obtained from 
function B,(T ) has significantly greater intensity than that 
obtained from B,(T ). Consideration of Fig. 3 leads to the 
same conc hision, the figure showing the curves for B,(T ) 
and B,(T ) constructed with different scales along the y 

axis, The ratio of the scales equal e2/z, and becomes huge 
for z= 5, Therefore, the integral of function By(T) will 

be much larger than the integral of B,(T), i.e., 


W (/)* = Wy, (/)* + Wa (/)* = Wi (f)’- (13) 


Since W,(0) = 0, W,(f)* = W,(f), and the spectrum of 
(7) becomes equal to 


W (/f) =W2(0) + Wi (/), 


f 
a where W,(f) is determined from Formula (12), 
The analytic expression for the spectrum 
W, (0) = 4\ B, (x) dt (14) 
0 
for5< zs 0, has the form (Cf. Appendix II); 
W, (0) = 29Af a (15) 


For telemetering systems, Fr/Af = 0.2. For fyjg,/Af =F)/Af <= 0.2 andn=1 
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Fig. 3, Correlation functions B,(T ) and B,(T ) at the discriminator 
output, 





exp |— 5 (“Ey] ~1. 


and in the range 0 = f = f,,9y. Formula (12) takes the form: 








_ 2n®Af / f 2 
Wi = =(a) (1+ ava. + ayaa: si 


For z2 5 








_ 2mtAf/ f\2 9.240 / 
Wi (a) =< (a) - pies 


The total noise energy spectrum at the discriminator output for fax /Af = F-/Af = 0.2 and 5 = z = 201s, 
according to (17) and (15), 


19.74 — 
W (/)\=— wee ty + 29Af =. (18) 


The noise power at the output of the low-frequency ideal filter 4g, with a receiver band-width of 0 to 


F 
W.= ( W (af = 6.58 22 (TE) + 29 ap? (54), (19) 


and the effective noise voltage is 


VW, = [6.5842(7t FY + 20ap =(H)] - (20) 














3. Optimal Frequency Deviation 





By definition, the magnitude of the mean square error at the output of the frequency receiver equals 


Oo == Dy a ’ (21) 
where U.,ay 1s the maximum dc voltage of the signal at the receiver output equal, according to Formula (1), to 


WD max» Where Wp max 18 the maximum frequer zy deviation. 


To determine the error, we express the quantity z in terms of the signal voltage and the specific noise 
voltage, Oo: 





:= (7) cer (22) 


Substituting the magnitudes of VW and U;j4x in (21), and replacing z by its expression in terms of 0 from 
(22, we get 


u2 "I, 
2 FP Fe os =% Af 
oo o.089(¢-) Aa + 0.26 (iE) ane 7 sol. (23) 


D max 


Formula (23) can be transformed by the introduction of the generalized parameters 


s 2 Af “Dma x 
= —, = —andy = 
P oV Ff : !D max . 








With these parameters, 


e* 


“a 
. 0.083 0.26 -* 
8 | sa ¥ Bying ° *) 53 te 


Consequently, the expression for the meansquare error 5 consists of two terms, Both terms decrease with 
increasing g, but behave differently with increasing y, i.e., with increasing frequency deviation, The first term 
decreases with increasing y, but the second grows. For large 8, the second term is much less than the first, and 
Formula (24) coincides (to within a factor) with the formula obtained for the mean square error from the theory 
of potential noise stability. 


It is possible to:determine from (24) the magnitude of the frequency deviation for which the error will be 
a minimum (considering n = const), We equate the derivative with respect to y to zero, thus obtaining the con- 
dition for determining Y opt’ 


p* 





‘. opt” (i pbs 1) = =e . (25) 
For z ze. 5, own 1, and Equation (25) will take the following form: 
Yopr! Yopt? 
o- Fope = 12 (28) 
Ben 











From the transcendental equation thus obtained, we can find Yopt AB» 1), 1.e., we can determine the 
dependence of the optimal frequency deviation on the signal voltage u,. the specific noise voltage 4, and the 
filter output band-width Fr, determining the system's rapidity of action. 


As is clear from the results obtained above, Formula (26) is correct for 





9S 2 on go Wand “t og < 0.2 
> “opt 27 opt” Ss A/opt Yopt” U.S. 


We see from (26) that 8 = g < 250. This condition determines the limit of applicability of Formula (26), 
Knowing Yopt» one can determine the least error possible with given values of g and n. Substituting (26) in (24), 
we obtain 





_ 70.083 , 0.33 4% 
omin= roa + PY opit | 7- (27) 


It is obvious from (27) that, for given 8 and y¥9_4, 5 min is minimum for Af/fp ,,,, = 1, i.e., when the 
maximum frequency deviation equals half the receiving band-width at the filter input. 


, a Figure 4 gives the dependence of yy, and 6,4, (for 
8 Yop” ©8 © min Jopr? n = 1)0ng, computed from Formulas (26) and (27). 
3 
a It follows from the curves that, for example, for g = 10, 
2 | 19 = 6.3, and for Fr = 0.5 cps, the deviation fp max" 


7 Yop 
> ae = 3.15 cps, but the error 6.4, = 0.56%; for g = 50, Yopt? = 
/ N oA = 100, and for F; = 0.5 cps, the deviation fp max = 50 cps 


q 
but the error 6 = 0.7 x10 
he log 8 —_ ™ 


5 Ag d It is clear from Fig. 4 that these relations are almost 














linear when plotted on a logarithmic scale, For 8= g =< 250 
(the range of practical variation of this parameter) these 
-2 \ dependences are very well approximated by the functions 








log (Yop) =—1.0 + 1.8 log 6, log (8 min)n=3 * 
= 2.8 -2.9 log B. (28) 


























= 
anit”! Here, 5 min is expressed in per cent. Dependences (28) are 
given in fine lines on Fig. 4. 


Figure 5 gives the functional dependence 6 = ( Yn=y’ 
for g = 15, 25, 40, all computed from Formula (24), It is clear 
Pp 25 50 75 0 100 from the Figure that, for a definite value of the frequency 
As y= 1D max deviation (y = Yopz), the mean square error at the frequency 
\/ Ff receiver output is a minimum, With increasing 8, 5 pin is 
SY 
N“ 








oe decreased, and the value of y for which 6 =6_,;,, becomes less 





sharply defined. 
4, Experimental Results 


| Experimental verification of the dependence of the error on 
the frequency deviation were carried out under laboratory condi- 
Fig. 5. tions, For this, use was made of a frequency receiver with a two- 
way discriminator for strongly limiting signal amplitude, The 
mean square error was measured by a special instrument developed at IATAN Laboratory 4. An unmodulated 2 
kilocycle signal was applied to the receiver input, as well as fluctuating noise from noise generator LIG-24-48, 


























Figure 6 gives the dependence of the error on the frequency deviation, as obtained experimentally (1) and 
by computation (2). 








It was assumed for this that U,/o = 36, Fy = 4 cps, 






































b% ; fp max = Af. 
a3 ; =. It is clear from Fig. 6 that there is good coincidence 
Sibi between the experimental and computed values for the 
a2 . A optimal frequency deviation, for the minimum error, and 
RL. Fa the character of the dependence of the error on the fre- 
gl tome quency deviation, 
0 BD HS 75 100 fp» ¢Pps SUMMARY 
Fig. 6 1, For a telemetering system there is an optimal 


value of frequency deviation for which a minimum error 
is obtained from the action of fluctuating noise. 


2. The optimal frequency deviation is determined by the value of the generalized transmission parameter, 


B = U,/o VFg For 8 = 8 = 250, the optimal frequency deviation (y,,, = f /F,) and the minimum error, 
Pp’ D max 
5 min» are determined from Formulas (28). 


The author wishes to express his thanks to G. A. Shastov for a series of good counsels, given to aid in the 
performance of this work. 











APPENDIX I 
Energy Spectrum W,(f) 
To find the spectrum W,(f) from (9), we simplify the expression for B,(T ) from (4). 
It is well-known that for large x 
- “1 > # (n+1)! 
E; (2) = 258 +rn ), nl@) << oy : a. 1) 


For large x and a fixed, but not too large, n, tp(x) will be small in comparison to the other terms in the 
square brackets, For x = 5, 


e 


i 2 6 
If E(z /p) is determined from Formula (I. 2), Expression (4) may be transformed to the form 


2 
B, (t) = a{® (1—2at%) + © (t—2x 2a) + Fr 1—3x2ae9} = 
Mai (I. 3) 

=a> 2 = 1) 4 —n 2at%), 


Figure 3,a gives the curves of B,(T)z/a for z = 5, 10 and 20 (for positive T), We determine W,(f) by sub- 


stituting B,(T ) in (9). Integrating twice by parts for each term in the sum, and substituting the value a = n(2Aef, 
we obtain 








win = LLIN tla), ela, 2 ate) 





z 2Va2° 3V52 
antAf ’ f 1y (a—iyp = =-+2(21) (I. 4) 
-=7 (ar) Aare ete 











APPENDIX II 





Energy Spectrum W (0) 


For finding the spectrum W,(0), we consider the graph of the function B,(T e/a , constructed from the 
exact Formula (5) (Fig. 3,b). To find W,(0) from Formula (14), it is necessary to find the area of the figure bounded 
by the graph of the function and by the coordinate axes, Simplifying Expression (5) for By(T ) over the entire range 
of variation of T (from 0 to oo) is difficult, since the area can be determined only by numerical integration, How- 
ever, for Var» 0 and Var — o (i.e., when it is difficult to obtain the necessary accuracy by numerical 
integration methods), an analytic expression for B,(T ) is readily found. 


We consider these cases, 
The case when at* << 1, Here, 





—pe at? 
= [oer twats. (il. 1) 


e=e * —4—ar?, 





If, in addition, a f*z << 1, then Formula (5) may be simplified, For small x 








E; (2) = C+ine+ 3) (C = 0.577). (Il. 2) 

n=i”™ 

For x << 1, E(x) *C + Inx + x. 
Then, 

B, (pat =) —26, (=? +) = —c—im 20%, (il. 3) 

2z 
2e~ “site, ~ 148 ~_ ¢-* (44 2). (Il. 4) 

~t—o i—po 


Substituting (II. 3) and (II. 4) in (5), we obtain 


By (t); = ae? {— 2az* (4 + 2) — (1 — 2at*s) (C 4- In 2at*2) \ . (I. 5) 


The case when af* >> 1, Then, p<< 1 and, for z = 5, it may be considered that 











—_ —z 
=P e>1, tee ~¥46 < mS 
and 
i—p z\ i—p Il. 6 
B (ize 5) < 284 ( e ‘). =s 


Substituting into (5) the value of , aha ‘ determined from Formula (I, 2), and discarding small quan- 
tities, we obtain, for z = 5, 





Paton ae te Lamas — Te]. au.) 


From Formula (Il. 5), one may determine the value of Va 1, (Fig. 3,c, curve for z = 5) for which the 
function BT); € "/a, becomes equal to zero. Setting By(T ), from (II. 5) equal to zero, we obtain, for z = 5, 
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(II, 8) 


eet 
For z= 5 
nd att =0,033 <1. at?z = 0.165 <1. (II. 9) 
> Consequently, the area of the region from VaT = 0 to V@Ty (Fig. 3,b, striped portion of the curve for z = 5) 


From the graphs constructed for z = 5, 10 and 20 (Fig. 3,b), it was found that 
Vat, = 1.55/¥r. 
On the basis of this, we may write 
ast 


W: (0) = Wa (0), + Ws (0)7=4| Balt) de +4 
0 


td 


expressions for W.(0), and We({0),,. Finally, we obtain 
p ) y 























may be found by a definite integration of Function (II. 5), Formula(II. 7) may be employed for finding the area 
over the whole range VaT > VaTy, It was. found, by a numerical integration to deter-aine the area of the central 
portion of the figure, that the area bounded by the portion of the graph of Function (5) from V&@T, to@ was 
numerically equal to the area of the figure bounded by the graph of Function (II. 7) from V&T, to o., 


(II, 10) 


(Il. 11) 


Integrating Function (II. 5) between 0 and 1T;, and Function (II, 7) between T, and oo, we easily obtained the 


W, (0), =3.74V¥q — (ll. 12) 
z 
and 
oaat oo Sae~? 
0 he a o — oe TPs 
since integration by parts of the first term in brackets gives 
) _ eeu 
ete OF at ee .14 
s, (te 3 aes | ra etic 
For §5<4:<20,E= aS pie =0.33to 0.39. Using (II. 10) and letting € = 0.36, we obtain 
Ws (0),, = 4.46V2 2 (Il. 15) 
Vz 
) 
Substituting (II. 12) and (II. 15) in (Il. 11), we obtain, for 5 = z = 20, 
an- 
mY hs e2 
W(0) =8.2Va == 29Af>-. (Il. 16) 
It is clear from Formulas (14) and (II. 16) that, for Af = const, the intensity of the noise energy spectrum 
1) —% 
for f = 0 is directly proportional toe /V¥Z. 
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A SELECTIVE LOW~FREQUENCY RC~AMPLIFIER AS A 
CONTROL SYSTEM ELEMENT 


Iu. G. Kochinev 
(Leningrad) 


The dynamic properties of a selective low-frequency RC-amplifier 
with a double T configuration are investigated. 


1. General Considerations 





Two modifications of selective low-frequency RC-amplifiers have strongly recommended themselves in 
use [1]: a circuit with a Wien bridge and a circuit with a double T configuration (Fig. 1, a and b), The tuning 
frequency, in both cases, is determined by the relationship w) = 1/RC. Both circuits are capable of guaranteeing 
high selectivity. 
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Fig, 1. 


In the actual realization of these circuits, difficulties arise in summing the input and feedback signals. 
More simply constructed solutions are obtained by the use of double T configurations, We turn now to such a 
circuit, and investigate its dynamic properties, 


The general questions of matching the feedback loop parameters with the amplifier have been clarified 
in sufficient detail in [2, 3], and will not be treated here, 


The transfer coefficient of the unloaded double T equals the ratio of the steady-state output and input 
vectors when sinusoidal signals are passed: 


3 w? 
@® @) 





(1) 


rp oem wo? — 4jaw, — w? 
In order to obtain high selectivity, it is necessary that the amplifier, without the feedback loop, have high 


gain and comparatively little nonlinearity and frequency distortion, all within the given frequency range, On the 
other hand, the input and output impedances of the feedback path must be matched with the amplifier parameters, 
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Figure 2 gives a schematic of a circuit for a selective 
amplifier satisfying the conditions imposed, This circuit may 
be considered a development of the circuit described in [4] 
for a double triode with cathode followers. 











The high gain in the amplifier we are considering is 
obtained by means of a reduced potential on the anode and 
screen grid of the pentode, and by the use ofa high ohmic 
load resistance (a vacuum tube functioning in the so-called 
“starvation” regimen), The cathode followers (T, and Ts) 
implement the impedance matching of the selective circuit 
and the amplifier cascade. 





























We now determine the transfer function of a selective 
RC~-amplifier, constructed according to the circuit of Fig. 2. 
To simplify the argument, we introduce the following simpli- 























‘, fications, which will allow us to go over to the circuit of 
Fig. 3: 
: atlas 1) the cathode followers have unit gain; 
r 2) the cathode connections of Ty and Ty are of negligible 
. influence; 
Yin? *elMout 3) we will conditionally assume that the double T does 
+ not load the amplifier cascade, Such an assumption is per- 


missible since, in actuality, the selective circuit is connected 


Fig. 3. to the low-impedance output of cathode follower Ts. 


We assume that the gain of the unshunted amplifier is a real quantity and does not depend on frequency, 
Such an assumption is excusable in practice, both for unwelcome harmonics, destined to be suppressed, and for 
side frequencies arising from carrier modulation. 


We might mention that there exist circuits capable of amplifying signals of these magnitudes without 
distortion, thus allowing amplifiers to be considered as linear elements. 


Keeping in mind the remarks already made, and using the notation of Fig. 3, we may write 


UG =Uin + & (w) BU out» (2) 
Vou gt ae ko U, (3) 


where g is the feedback coefficient and-ky is the open-loop amplifier gain, 
Possible cases are g < 1 and g > 1; the latter case indicates that there is an amplifier in the feedback path, 
From (2) and (3) we obtain the expression for the gain of the selective amplifier: 


& (w) — You _ he (4) 
= Dig Tao phe 





2. Transfer Function of the Selective Amplifier 





We will call the amplifier transfer function the ratio of the’ Laplace transforms expressing the envelope at 
the amplifier output and input for zero initial conditions: 


FE out (P) 


: (5) 
Ejn(P) 


k™ (p) = 
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The carrier transfer function will be the name given to 


_ Yeu 


Aer (8) 


With these definitions, (5) is permissible only in case the carrier is not subject to phase shift and the form 
of the rounding at the output is independent of the initial carrier phase, 


The concept of the complex transfer function may also be applied to the selective amplifier: 


; E 
k (Q) = roa (1) 


corresponding to the steady state, and characterizing the amplitude and phase variation at the output for sinusoidal 
rounding at the input, 


In [5] is to be found an approximate expression for the complex transfer function of a selective amplifier 
for small values of the modulating frequency Q: 





k 
&(Q) = Thor ie 
where 
, kod 4 
> Deny’ ) 


Based on (8), the selective amplifier is equivalent to a first-order inertial link with time constant T° and 
gain kg. The envelope formed at the output of such a link by passing a unit sinusoidal signal is described by the 
well-known equation 


t 
Eou(t)=k(1—e *). (10) 


It should be mentioned that Expressions (8)-(10), despite their apparent obviousness, give an insufficiently 
accurate picture of the amplifier behavior during the transient response and distort the physical picture of the 
phenomena for established unit sinusoidal impulses, To corroborate this assertion, we consider the passage of 
the following signal by the selective amplifier: 


U in(t) = sin (cot + 90). (11) 


The Laplacé transformation of (11) is well-known [6}: 


$ SIN Po + Wo COS Po 











Uj) = P+ (12) 
The transfer function (carrier) for the double T and k (8) are obtained from (1) and (4): 
: Atel 
Aiaias s* + dees + wo? : (28) 
. a s* + 4a s + wo? 
k (s)= TH hs : 


s? + ts s+ w? (14) 
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Then, in accordance with (6), we obtain 


—ky sin g,+s*(@,cos9,-+ 4e,sin ¢,)+ s(w)sin 9, + 4a) cos ,) + w)cos p, 














U—, (s) = on 
out 1 + ko3 2) 2).[ «2 4wo 2 
« +0p-(s +77 epi ts ) (15) 
iow Sia $8in Go + 0008 o , ; koB_ __ 8 SiN Po + COS Po 
s* + w? 1 + ko O4 sta} 
The source, in the time domain, of the first term of (15) is the steady-state solution: 
Uout st (t) = — ko SiN (coot + 9). (16) 


To determine the source, in the time domain, of the second term, we employ an approximate method 
germane to systems with little damping [6]. 


For the case under consideration, in the expression 


$ SiN Mo + Wo COS Po __. 5 SiD Go + Wo COS Fo 


O+ esto} s* + 8a,s + wh (17) 





as the “small parameter” of the problem we may take the quantity 6 = 4/1 + kpg), characterizing, as is easily 
verified [4], the “damping” of the selective amplifier. 


The roots of the denominator of (17), for 6 = 0, assume the values s} 2 = Jw. To find the roots lying in the 
neighborhood of jw», we set 8 = q + jw» (q is a small quantity, of the order of magnitude of 5). 


Then, (q + iuyy + Sw» (q + juy) + us = 0. Hence, dropping terms of the onder of magnitude of & , we find 
that jw (2q + wy 5) = 0, or q = —w,5/2. 


Thus, improved values of the roots will be 


4, =0,(-—+ +i), ‘, = o,(—4-—i) for + <li 


Now, we decompose (17) into partial fractions: 


$ 81D Po + 9 COS Po _ a; + a2 


s* + 805+ oF s—s, s—8, 





and, finding the inverse transformations in the time domain, we get 


a; Qs ¥ *,! -> 
oa r+ ‘| =ae +ae =e Sin (aot + 9p). 


s—8, s—5, 





Thus the source, in the time domain, of U>,;(s) takes the following form: 





Vout (t) = — kysin (oot + m)(1— He *), (18) 
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It is clear from a comparison of (18) and (11) that only the envelope which is unconnected with the initial 
phase of the carrier, is subject to variation, It is possible to leave out of consideration the phase inversion of 
the carrier, which is a property of all amplifiers with nonreactive anode loading. 


The “independence” of the transmission of the carrier and the envelope allows the concept of the transfer 


function to be applied to the latter (Cf. (5)J. Its values are found by applying the Laplace transform to the expres- 
sions for the envelope at the input and the output: 








t 
Ein(t)=1, Four (t) = hy (1— -Bh ae *), (19) 
tT 
ae 1 1+Pitka 
Ein) =» Foul) = he (> — 79g 5ae) = he a 


p 1+kBp+t 0 pa+pt) ° (20) 


We obtain the results 


Fou (p) , ‘t+PizEKs 








(21) 


The qualitative differences between Expressions (21) and (8) are quite evident, The establishment of an 
envelope (21) with unit sinusoidal impulse is described by Equation (19), where the final value of the output 
signal, corresponding to the moment of switching, equals 


E out (to) = ko (4 > ka) - 
ko (22) 
T+ kS° 


At the same time (21) and (8) should not be considered as in contradiction to one another, Their identity 
for ky 8>> lis obvious, Expression (21) should be considered as the more accurate and correct reflection of the 
physical side of the process, 





Figure 4 shows the oscillograms of the process occur- 
ring when a sinusoidal voltage was applied to the input of 


$V" a selective amplifier built according to the circuit of Fig. 


2, using tubes 6N2P (tubes T, and Ts) and 6Zh1P (tube T3). 


U] 
out 
ee OV The values of resistance and capacitance used in the cir- 
*o 


cuit were as follows: R = 680 kilohms, Ry = 1 megohm, 


Us Rg = 10 kilohms, Rg = 1,6 megohms, Rg = 180 kilohms, 
Sf J IVW/VVVV Rg = Rg = 2.2 megohms, Ry = 20 kilohms, Rg = 270 kilohms, 
b ab C = 5000 pyf, Cy = Cy = 0.1 yf. 
B-Q'5, 6, +180 The upper curve corresponds to the input signal, the 


lower curve to the output signal, The form of the oscillo- 
grams, read off with two values of feedback (B, = 1.0 and 
B}, = 9.15), is in good agreement with (21). 


In conclusion, a remark of a practical character should be made, It is more correct to consider Ko, not 
as the open-loop amplifier gain, but as the gain at the tuning frequency. For large values of gain, these two 
quantities can differ considerably. 


Fig. 4. 
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PHASE DETECTOR FOR MULTIPLE FREQUENCIES 


R. la. Berkman 


(Lvov) 


In a host of circuits, the output quantity is a voltage whose frequency is an even multiple of the frequency 
of the exciting voltage. The most widely used circuits of that type today are magnetic amplifiers and magnetic 
probes of the "second harmonic” type, wherein the output quantity is a voltage of doubled frequency. 


The rectification of second harmonics by means of phase detectors of the usual type involves complicating 
the circuit, since a frequency doubler is necessary. Differential amplitude voltmeters, frequently employed for 
the same purpose, have poor zero stability, narrow limits of linearity and provide low output power, 


A new circuit for a phase detector, due to M. A. Rozenblat [1], is based on the use of nonlinear carborundum 
resistors with symmetric voltampere characteristics, This circuit is simple and possesses high zero stability, It 
is not free, however, of certain drawbacks, connected with the facts that carborundum resistors work reliably only 
with high controlling voltages and pass only comparatively small currents, The controlling power thereby ob- 
tained is considerable, Therefore, in those cases when devices of high efficiency are called for and significant 
rectification of current is required, the circuits mentioned are hardly suitable. 


It should also be mentioned that carborundum resistors of the necessary dimensions are not manufactured 
industrially at present. 


It is stated in [2] that instead of carborundum resistors, one may use reverse-connected semiconducting 
rectifiers, Such a solution to the problem is impossible to accept as satisfactory since, in the first place, the zero 
stability of such a device is extremely low and, in the second place, the curvilinear portion of the rectifying 
characteristic is very small, which introduces difficulties in the computation of the circuit (a large number 
of series-connected rectifiers is required). 


We have developed a new circuit for a phase detector for 
multiple frequencies, similar to the type of detector with non- 
linear symmetric resistances (Fig. 1). It is based on the use of 
nonlinear elements with symmetric voltampere characteristics, 
taking the form ofa bridge circuit with four rectifiers, so con~ 
nected as to cut off all voltage (in Fig. 1, the bridge is enclosed in 
dotted lines), With this connection, the nonlinear element pos- 
sesses a great deal of symmetry, since opening of the link occurs 
for one and the same value, u > Up, independently of the direction 
of the applied voltage (Fig. 2). The influence of other sources 
of asymmetry — variation in the direct resistance of the rectifiers — 
can be decreased by connecting additional resistance in the cir- 
cuit, 


The phase detector built on the circuit of Fig. 1 may be designed easily for variations, within wide limits, 
of the magnitude of the putput resistance, sensitivity and linearity with respect to the measured voltage, by means 
of a limited number of rectifiers, by varying the cutting-off voltage and the auxiliary resistance, It is essential 
for this that, if auxiliary resistances are introduced, the power in the rectifying circuit exceed slightly the power 
dissipated by the rectifiers, 








Fig. 1, 
































a = cot™*(2 r,) 
7 
4 +4) r 0 5 -- i] F 
Myr dug rr. * Yy 5 
Fig. 2. 
“Y, lee Sa e eee, an _ 
a Fig. 3. 

4, sige , 

- ~ Zwt Figure 3 clarifies the principles according to which 
” our circuit works, During the intervals of time Ty to T; 
it and T, to Tg , the circuit is cut off, since u< Uy. During 


the intervals T, to T; and T; to T,, the circuit conducts 
and, for a definite current i, the principle of superposi- 
tion may be used, The average value of the commutating 
current during the times (T; to Tg) + (Tg to T,) equals zero, 
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Fig. 4. 








Thus, the circuit functions as an interruptor with 
doubled frequency, where the “closed” time is controlled 
by the magnitude of the closing voltage U, (for a given u,). 


For the analysis, we assume a piecewise linear voltampere characteristic approximating that of the rectifier. 
Detecto; behavior is investigated during the steady state. We may write the initial relationships: 


commutating voltage 
Ue = Ucmsin oat; 


measured voltage 


— TU man) sin (2net — Pan) + U nian) sin [(2n + 1) wt — Panta) 


n=1 


closing voltage 
U, = En (1 + 8) + 2U or; 


total resistance to ac in the circuit 
r=2rp+7riztrictry 


We have used the notation: Up, is the initial voltage across the rectifier and r, is the direct resistance of 
the rectifier, 


We choose a nominal value for the closing voltage Ej, starting from the condition that the circuit be in 
the open and closed states for equal intervals of time. This corresponds to 


4 U 
Ey= FF om UY, U,= 7 (1) 


or’ 


ony 


The quantity 6, characterizes the deviation of Ey from its nominal value;5,2 0,65. << 1, 
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The rectified current in load resistor R is determined from the equation 


Tt. ™%; 


k= + (\ idt + | iat) = 
= 5 [\ we +u,—U,—h,R)dt + ( (we +u,+U,—1,R) dt] (2) 


We determine the moments of time 1 for small measured voltages, neglecting the quantities u, and I, R,. 
Then 


U,, Sin or, = +0, (1 +8,), simmer, = + vy tt +8). 


Expanding arc sin wT; in a Taylor series, and limiting ourselves to the first two terms, we get 


bey 3 
or, = 7 + bs, ig ae ots = + +, on, = 2 —8,. (3) 
Since Ts = %] + t/w and Tg= T+ tho, 
\ (ue —U,) at + \(we+ U,) dt =0, 
\ D Umcon4aSi ((2n + 1) ot — 9, dt + 
a — (4) 
+} Sy Fonconss8it ((2n + 1) ot —9,,4,)dt = 0, 
™% n=1 


i.e., the odd harmonics of the controlling voltage are not rectified, For even harmonics we obtain 








tr 2 


L=- un (+ — 28,) — 3-72" cos (2nust — ,)| : 
" 


n=-1 


or, after transformation, 


. > a (sin Pon, 00S 2n8, sin > + sin 9,,, sin 2n8, sin =“ is *) 
= , 
*[r++(1-4-%)} (5) 


It is clear from the formulas we have obtained that, for the nominal value of closing voltage (5. = 0), the 
only harmonics detected are those with ordinal numbers in the series 2n = 2 + 4m (m = 1, 2,2.) 





For the second harmonic 





R sin ?2- (6) 


In the unloaded state (for R= o) 
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At maximum power in the load, determined by the condition that 6 P,/ OR = 0, we have 
Ropt = 2r. (8) 


Clearly, an optimal value of r does not exist, r,,,—» 0. The minimum value of r is limited by the 
- Pp -~_ 
admissible power required of the controlling voltage source, 





For use of the detector in a noncompensated circuit, an important part is played by the degree of nonlinearity 
of its voltampere characteristic. To determine the degree of nonlinearity analytically is, in the general case, 
quite difficult, since it is necessary to deal, in the expansion of arc sin wT, with quantities of greater than third 
order smallness (Cf, [3]). Linearity appears to be best for short-circuited states and for loads which are not shunted 
by capacitance. In these cases, as is clear from the physical picture of detector behavior, linearity is maintained 
down to the value of U,,,, for which the commutation angle (3) is disturbed. As is clear from Fig. 4, this cor- 
responds to the condition 


U xm<Uc m (1 — 73) or Ux<0,3u, for = 90°. (9) 


With the introduction of load resistance, the linearity deteriorates somewhat, as compared to the unloaded 
regimen, 


For the practical execution of the circuit of Fig. 1, it is convenient to replace the source of closing voltage, 
Ey, by a resistance shunted by a sufficiently large capacitance. The value of this resistance is determined from the 
condition that the voltage drop across it from the rectified commutating current will equal the required value of 
closing voltage: 


2° 
Ey Ro =F; \ (we—U,) dt, 


1 


fromm which 





2 
Ry = > 7.37. (10) 


A circuit for the detector employing capacitance-ohmic closing is given in Fig. 5,a. Resistor R., serves 
to establish zero. 
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Fig. 5. 


It is undesirable, in many cases, that commutating current pass through the source of rectified voltage. 
To avoid this condition, it is possible to use the symmetric detector circuit given in Fig.5,b. The sensitivity 
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threshold of this circuit is determined by the zero drift during 
— = the time of rectifier aging, variation of commutating voltage, 
ss L’ etc. Zero instability was experimentally determined for 





detectors of germanium diodes of type DG-Ts over a period of 

one day and, for variation of + 20% of the commutating voltage, 

‘\ was between 5 and 20 millivolts of equivalent rectified voltage, 
depending on the type of rectifier and on the matching of 

100} Q5 \ characteristics, 


Be inet Thus, with respect to sensitivity threshold, our detector 
does not differ from the usual circuit using semiconductor 
rectifiers, Its use seems expedient in those cases when the 
j “; Volts attainment of a minimum sensitivity threshold is not essential, 
° ” 7 od %~VO\S for example, in devices with second harmonic preamplifiers, 
Fig. 7. in intermediate and output cascades of magnetic amplifiers, 
etc, 
































In detectors of the type we have been considering, there is also rectification if the frequency of the com- 
mutating voltage exceeds that of the measured voltage, For an example, we consider the behavior of the circuit 
of Fig. 1 with the conditions 


Uc=Ucmecosat, Ux = Usmsin >t, UjnUe = Se. 


11 
VE (11) 


It is obvious from Fig, 6 that, to determine the magnitude of the rectified current, there should be carried 
out an integration over two periods of controlling voltage: 


Ts 


h = ge [\ (uc + ue—Uy— A R)dt +\ (te + te + 0p — MR) dt + 


lait vt mn (12) 
+ | (we + tz — Uy — HR) dt + (ue + tx + Uy—¥ R) at |. 


We compute the values of T by means of Taylor series expansions of the functions Vy and arc sin y 
obtained from the previous solution, and retaining quantities up to the second order of smallness. 


The momentsof time 14, T,, Ts, ™ are determined from the conditions 

































Ux + Uc = U4, tz =Usm WV 4(4 — cos wx), 


Uc m (cosux — =) = —U xm V +a — cos wt), 


U2 nV EAU nV om V2—V2— +03, 
Qu, an ' 





cos @t = 





where we obtain 


m= —F4+8V2—RV248, oy = 2-4 V2—2, 248, 








(13) 
on= $43 V2+RV2—8,  oy= F— 2, 7244/28, 
where 
g, = V2— V2 Vem 3, — 1 Vem 
Oks Mam “Th. 
Analogously for Ts, Tg, T;, T,, we find 
ory + mVAFRVI+8, w= 3,243,248, | 
(14) 
where 
Ries V 24+V2 Yom 
2 cm 
Transforming Equation (12), we finally obtain 
u2 


: aa ress > 


It is interesting to note that in detectors with such circuits, as opposed to detectors of other types (including 
detectors of carborundum resistors), the magnitude of the rectified currentis inversely proportional to the com-~- 
mutating voltage. This might be exploited in a host of automatic devices, 


Experimental curves for ly, = f;(u,) (curve b) and I, = f(u.) (curve a) are given in Fig. 7. Curve a is given 
for u, = 2 volts and r = 1 kilohm, curve b is given for u. = 30 volts and r = 1 kilohm., 
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CONCERNING SOME PROPERTIES OF FERROMAGNETIC CLUTCHES# 


P. N. Kopai-Gora 


(Moscow) 


Ferromagnetic clutches are compared with electric motors and 
friction and hydraulic clutches in terms of their starting time, electro- 
magnetic time constant, control power, and theratio of the output 
moment to the moment of inertia. The said parameters of the ferro- 
magnetic clutches are shown to be better than those of the electric 
motors and friction clutches. 


Lately, the transmission of torque and its control is being accomplished by means of ferromagnetic 
powder clutches. The working space in these clutches is filled with powdered carbonyl iron together with a 
dry or liquid filler. Such a mixture, when acted on by a magnetic field, is able to couple the surfaces which 
enclose the working spaces, thereby permitting the transmission of torque between the shafts, 
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Fig. 1. Dependence of the transmitted torque on the 
excitation current and speed of rotation. 


The dependence of the transmitted torque on the excitation (a) and angular frequency (b) for ferromagnetic 
clutches is shown in Fig. 1. The diagrams show that the clutch characteristics are very favorable for use in auto- 
matic control systems, since the moment is a linear function of the excitation current over a very wide range, 
and is independent of the speed of rotation. 


The following qualitative comparison of such clutches with electric motors and friction and hydraulic 
clutches is carried out to establish more fully the possibilities of utilizing ferromagnetic clutches as actuating 
devices in automatic control systems, The comparison of the indicated devices is carried out on the basis of 





their starting time, control power, and the relationship = a (where M is the output moment, J is the 


M 

2 vy 
moment of inertia). 

2 There is no standard name for such clutches at present. They are called powder clutches, clutches with a 
magneto-liquid or ferromagnetic filler, etc. 


® Motors are usually evaluated in terms of the relationship M/J, but such an evaluation does not always prove 
to be satisfactory since the system's properties also depend on the transmission ratio between the motor and the 
load, 





The value of « for a motor of the SL (a) and PN (b) type, depending on the power, rpm, etc., varies 
between a few units and several tens, as is apparent in Fig. 2 (the numbers adjacent to the points indicate the 
type of electric motor). 


The comparison is conducted with the following assumptions: 


a) the mean induction in the gaps  —— of the samples of the ferromagnetic clutches and motors 
being compared are equal; 


b) the dimensions of the rotors of the electric motors and ferromagnetic clutches, as well as their moments 
of inertia, are respectively equal; the dimensions and moments of inertia of the friction discs and the ferro- 
magnetic clutches are also equal; 


c) the design coefficient of the polar arc («') for dc electric motors is equal to the ferromagnetic clutch's 
coefficient of utilization of the active surface (k). 


Such a comparison permits, in the design of automatic control systems, the direct choice of the actuating 
mechanism from the torque, and also permits an indirect comparison of mechanisms having equal power 
ratings if the angular speed is known* The mean tangential force per unit area for ferromagnetic clutches, over 
the linear portion of the characteristic (for 500 < B,,¢,, < 10,000 gauss), and for the usual relationships between 
the filler and iron in the work mixture, can be approximated by the formula 


T = Bmean Dk x 10“ kg/cm*, (1) 


where n is the number of pairs of active surfaces or the number of work layers of the mixture. 
In the case of dc electric motors the mean force per unit area p can be approximated by 
P = 1.02 Brean Aa'+ 107" kg/cm’, (2) 
where A is the linear load on the electric motor. 


Since, according to the above assumption, the mean inductions and moments of inertia are equal and, 
moreover, a’ =k, the following relationship is valid: 
ae n-10~4 n 


T 
z=-s e— 10°, 3 
Gir 2.) LOA” © <A , @) 





Since the magnitude a is independent of the transmission ratio and characterizes the ability of the drive 
unit to impart an acceleration to the load, then, in those cases when the transmission ratio can be chosen, it is 
more convenient to evaluate the motor in terms of the magnitude a = MAJ . 


fe ct for high speed control systems, the acceleration of the load shaft is, ignoring friction, 
s Tari an he where i is the transmission ratio between the motor and load, Mg, is the driving torque at the 


motor’s shaft, Jgy and J; are the moments of inertia of the motor and load,respectively. Differentiating this 
expression with respect to i and equating it to zero yields the optimum transmission ratio which provides maxi- 
mum acceleration. The expression for the maximum acceleration of the load shaft will be 
€ = isa a =a : 
ae dr l “ vy, 


Thus the bigger a is, the greater is the attainable acceleration of the system's output shaft. 

The mean value of the magnetic induction in the gap is understood to be that value which, when multiplied 
by the gage length of the polar arc, is numerically equal to the integral of the actual induction curve taken over 
the whole length of the polar arc. 

4 The comparison of ferromagnetic clutches with other types of actuating devices having the same power 
rating is complicated by the fact that, at the present time, the optimum geometrical dimensions of certain 


types of mechanisms have not been developed (for example, hydraulic mechanisms, ferromagnetic and friction 
clutches), Moreover, the power developed by the clutches can vary over a considerable range, depending on 
the nature of the cooling system. 
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Fig. 2. Dependence of a on the power: a) for series SL electric 
motors, b) for electric motors of the PN series. The numbers 
adjacent to the points denote the type of motor. SL~569 is de- 
signated by 569. The signs used for the points in Fig. 2,a cor- 
respond to the rpm's of the motors: @) 850, O) 2400-3800, 

X) 3000-3800, 0) 4000-4500 and @) 5200 rpm. In Fig. 2,b the full 
curves join points corresponding to one type of electric motor. 
The point designation corresponds to the motors’ rpm's; @) 1000, 
X) 1450-1500, O) 2000 and @) 2800-2870 rpm. 


where a, and a. .545, are the M/VJ expression for the ferromagnetic clutch and electric motor, respectively. 


The linear load on conventional electric motors, depending on the power and other factors, is in the region 
of 100 < A< 450 amp/cm, where the upper values of A apply to high-power electric motors. 


Then, for n = 1 (for the clutch shown in Fig. 3a) a, = (2.2-10) amotor- 


Drum-type clutches with hollow cylindrical rotors and having two pairs of active surfaces (Fig. 3b), are 
compared with direct current electric motors having hollow rotors. 


The moment of inertia, referred to 1 cm? of the cylindrical surface of the rotor, is expressed, for both 
devices, by the same formula: 
y 4D 
fa (4) 
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Fig. 3, Ferromagnetic clutches: a) usual type with one pair of active surfaces, 
b) hollow-rotor type, having two pairs of active surfaces. 





Here D is the rotor's diameter, y is the density of the 
material, g is the acceleration due to gravity. 


The linear load on an electric motor having a hollow rotor 
can be expressed in terms of the current density j , rotor thick- 
ness h, and the conducting efficiency of the rotor body, n 


=D aD = *U- (5) 





Then the force per unit area for such a motor will be 
P © Brean hnj*107" kg/cm*. (6) 


Here | is the total rotor current and S is the cross section 
Fig. 4. Cross section of magnetic circuit sda dikes 
for a high speed ferromagnetic clutch. : 





Using expressions (1), (4), and (6) we can write 


ae _ (Bmean kn), 10° als (hy motor ; (7) 
&motor (Bmean 12 j)motor (hy D*), 











Taking into account the equality in geometric dimensions and mean induction, and taking j = 500 amp/cm?, 
k =n = 0.7, n = 2 and ¥¢ = Y motor WE obtain 


4 2 x 10° ri ico (8) 
= to =———— r 
500 Dmotor — Bnotor eee 


The performed calculations are approximate but, nevertheless, they offer convincing proof of the pre- 
eminence of ferromagnetic clutches in comparison with direct current electric motors. The indicated relation- 
ships remain valid for alternating current electric motors. 


Let us now turn to a comparison of ferromagnetic and friction clutches, At the same time we will assume 
that disc clutches are being compared. The unit tangential force q for friction clutches is, as is well known, 
expressed by the formula 


q= 0.4 _ fe ny x 1077 kg/cm*, (9) 
Here f. is the coefficient of friction, n, is the number of pairs of active surfaces. 








Taking into account that the moments of inertia for the discs in both devices are the same, we have 





-4 
a, _7 _ _ (@meankM)e 10 . kn x 10° (10) 
afc, 1 MBean fr MFC, 10°" 0.4 Binean/ fr ™ 


From here, if we assume Bmean = 8000 gauss, fs, = 0.2, k = 0.7, n = ny, we will obtain 


a 
c 0.4 x 8000 x 0.2 





a OF oc, = 8.7 bat 7 (11) 


Thus, on the basis of the value of a, ferromagnetic clutches are superior to electric motors and friction 
clutches, 


Example: Let us compare a piston-type hydraulic actuating mechanism with a ferromagnetic clutch, Let the 
actuator’s piston-travel be L = 15 cm, diameter of the piston 3.5 cm, pipe-line diameter d = 1 cm, pipe length 
1 = 150 cm, oil pressure py = 10kg/cm* . For these parameters the pressure of the oil in the cylinder on 1 cm? 
of the piston’s surface gy = 15 x 10°? kg/ cm and the pressure in the pipe-line on the cylinder surface is 

1.82 kg/cm*. The total oil pressure on the surface of the cylinder is gy = 1.835 kg/cm* . The value of a 
for such an actuating device is 


1 


r 
a al — (f) = 23 kg-m¥*. sec“, same 


Here r is the radius of the lever or gear on which the rod acts, g is the acceleration due to gravity. 


For a two-layer ferromagnetic clutch with a hollow rotor (n = 2; see Fig. 3,b) we can take T to be 
0.7 kg/cm* , k = 0.7, and the pressure of the rotor on its surface, for a wall thickness h = 0.25 cm, is equal 
to gg = 0,002 kg/cm*, Then 


tnkr g ~1 
rie ay V £=s% kg - m2 « sec“3, (1b) 


F V & r2 
€ 
The calculations show that, in this particular case, the ferromagnetic clutch is, on the basis of the value of 
a, 30 times better than the hydraulic clutch. However, it should be noted that for other types of hydraulic 


actuating devices (for example, the rotating type) the magnitude gs can be smaller by an order of one or two, 
while the working pressure can be chosen to be 100 atmos and greater, 


In connection with this, a for such actuating devices can take on values greater than those for ferro~ 
magnetic clutches, For the comparison of concrete types of ferromagnetic clutches and hydraulic actuating 
devices in terms of a one can make use of the formula obtained from dividing (11b) by (11a): 


Let us now consider the starting times of various devices, The starting time t, can be approximated by 
the formula® 
“a 
Jdo Jn 
ts= \ 


= ’ 


: 6Mn = =BMn (13) 


® The example is taken from [2]. 


¢ Equation (13) is a very approximate one, but its use in this case is justified to the extent that this comparison 
of the various actuating devices is a qualitative one. 
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where g is the multiplication factor of the starting moment, which here is assumed to be equal to 1,5 for all 
the devices, M, and wy are the nominal values of the moment and the angular frequency, respectively. 





If in (13) we express M, in terms of a, we then obtain 
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In accordance with the above assumptions the multiplication factors of the starting moments and the 
moments of inertia of the devices being compared are respectively equal. If we further assume that the 
nominal speeds are equal, then 


Ao 
» ‘sy = ay &, - (15) 


Thus the starting times of the actuating devices being compared are inversely proportional to a. For ferro- 
magnetic clutches the starting time turns out to be from 2 to 20 times smaller than for electric motors and 
several times smaller than for friction clutches, The starting times of hydraulic actuating devices can be, 
depending on their construction, smaller (for high pressure, rotating actuating devices) as well as larger (for low 
pressure, piston mechanisins with a long piston travel) than those of ferromagnetic clutches, 


Let us compare the electromagnetic time constant of the ferromagnetic clutches to the time constants 
of other mechanisms. 


The magnetic circuit of a ferromagnetic clutch can be achieved by means of several configurations, One 
construction of a magnetic circuit for a fast-acting ferromagnetic clutch, containing an iron mixture, is shown 
in Fig. 4, The layout of the magnetic circuit of such a clutch is very similar to a direct current electric motor 
with the sole difference that the air gaps between adjacent poles, which exist in electric motors to reduce 
leakage currents, are absent in the clutch. This permits better utilization of the stator’s periphery and the choice 
of a large number of poles. We will show that the electromagnetic time constants of a single-layer, ferro- 
magnetic clutch and an electric motor's field coil associated with one pair of poles, are equal if the dimensions, 
gaps, and inductions are equal. At the same time we will ignore the leakage currents and ampere-turns, which 
cause a magnetic flux to be set up in the steel, in comparison to the ampere-turns necessary to set up a magnetic 
flux in the gap. 


Let Z = #D/p be the pole width, where D is the rotor's diameter and p is the number of poles. 
Then the inductance of the field coil on one pair of poles will be 


(16) 


Here I is the current in the field coil, W is the number of ampere turns, § = ZI is the pole area (2 is the 
pole length). 


It is known that the magnetic induction, ampere-turns, gap size 6, and magnetic permeability y are 
related through the expression 


[Wy 
Bmean= 95 ° (17) 
The coil resistance is approximated by the formula 
l 2w 
Rap lt22W = 


% 


















where q, is the wire gage, p is the specific resistance of the wire. 


Taking Equations (16) to (18) into consideration, the electromagnetic time constant of the coil will be 


_ BmeanS%s _ _u WSq, . 
Qp(l+Z) 46p (1 + Z) (19) 





L 
Tem =R 


In Equations (16) to (19) all quantities, except » and W, in accordance with the assumptions made 
earlier, are equal for both devices, 


As a result of the equality of the magnetic inductions, and taking into account the fact that the magnetic 
permeability of the mixture y1;, 1s 5 to 8 times greater than the permeability of air y,, we can write 


TWelt mix = Wma (20) 


If Jm = Jc, then —, +m = 5-8 and, therefore, the product 4 W remains constant in (18) for both devices. 
a c 
In this way, the time constants of the ferromagnetic clutch's control winding and the electric motor's 
field coil are equal when their dimensions, gap sizes, and magnetic inductions are equal. 


The control power for ferromagnetic clutches is 5-8 times smaller than the excitation power for a direct 
current electric motor since yp, = 5-8 1.4. Taking into account the fact that the excitation power comprises 
only 3-10% of the motor's total power, it should be noted that the power required to control a ferromagnetic 
clutch represents only 0.6-2% of the transmitted power. 


The small controlling power required for ferromagnetic clutches permits the use of quick-response, inertia~ 
less, control units (e.g., electronic) to secure the forced operating mode of the actuating mechanism and to raise 
the quality of the whole system. 


The above comparisons of the various devices specified equality in their dimensions, i.e., the equality of 
the active surface. Obviously, the ferromagnetic clutch can, with the help of great pressure, develop from 2 
to 10 times more power than the electric motor. For the same power, the active area DB in the ferromagnetic 
clutch can be chosen to be 2 to 10 times smaller, which results, approximately, in an equal reduction of the 
time constant and clutch control power. 


It should be noted as well that ferromagnetic clutches can be made with a large number of poles, which 
also results in a significant reduction in the time constant. 


The considerations presented here show that the time constant of ferromagnetic clutches will be consider- 
ably smaller than the time constant of an electric motor of the same power rating. This is confirmed by experi- 
mental results, Thus for a clutch with M, = 10 kg-m the time constant turned out to be equal to 20 msec, 


The above considerations are also quite applicable to the comparison of ferromagnetic clutches with 
friction types, which leads to the following results: for clutches having equal dimensions, gaps, and magnetic 
inductions, the time constants of the coils are equal, while the control power for the ferromagnetic clutches is 
5 to 8 times smaller than for the friction ones, 


The dynamic properties of ferromagnetic clutches are determined to a significant extent by their cons- 
truction, Since in such clutches the torque is determined by the magnetic flux, designs aimed at decreasing the 
response time must take into account means of ensuring a quick build-up of the flux, In particular, to 
minimize eddy currents, the magnetic circuit should be made of a blended iron having a low magnetic resist~ 
ance and a high resistance to current; the elimination of all short-circuited turns must be ensured, Thus, for 
example, the clutch’s rotor should be a split type. Clutches with a liquid filler have a shorter response time 
since the motion of iron particles in a liquid medium is much more free than in a dry powder, To decrease the 
electromagnetic time constant it is desirable to design the clutch to have many poles, which permits a 
reduction in the volume of iron involved in a pair of pole pieces, 


However, execution of the indicated recommendations does not always guarantee the proper fast response 
of the clutches. To confirm this let us turn to some experimental results, Experiments were carried out on 4 








Fig. 5. Split rotor of a two-layer ferromagnetic clutch. 


cylindrical, two-layer, ferromagnetic clutch with rotor 
dimension D = 150 mm and width of active area 65 mm, 
designed for a moment M = 10kg<m, The rotor in the 
moment clutch was a split type, as shown in Fig. 5, and the magnetic 
current circuit was ablended one. As experiments showed, the 
magnetic flux practically coincided with the electric 
Aig msec current, One could suppose that in such a construction the 
f* cycles moment must also coincide with the current, but experi- 
mental results have shown that the moment lags behind 
the current quite appreciably (see the oscillogram in Fig. 6). 












Fig. 6. T is the build-up time of the moment The lag of the moment build-up behind the magnetic 
in 0 cinch itd uneiaied dom. flux is explained by the motion of the iron particles, under 

the action of centrifugal forces, through the rotor slots, 

from the inside gap to the outside one. When voltage is 
applied to the control windings the iron particles return, under the influence of the magnetic flux, to the inside 
gap. This process occurs in a certain finite length of time. Moreover, the steady state value of the moment sets 
in after all the particles have migrated. This last is confirmed by experimental results; after the slots were 
sealed the build-up of the moment in the clutch practically coincided with the current, as is apparent from the 
oscillogram in Fig. 7. 
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Fig. 7. T is the build-up time of the moment in a clutch with sealed slots. 


The clutch was a four-pole one, with equivalent gap 26 = 0.6 cm, Bmean = 7000 gauss and k = 0.5. The 
control power did not exceed 50 watts which permitted the use of an electronic control unit and guaranteed triple 


or quadruple forcing. The electromagnetic time constant of such a clutch turned out to be 15-20 msec which 
is less than for electric motors of the same power. 
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Ferromagnetic clutches do have some shortcomings. 


Under the condition of continuous slipping the clutches generate considerably more heat per unit area than 
electric motors; as a result considerable effort must be put into its elimination, 


Under high speeds of rotation centrifugal forces disrupt the normal operation of the clutch and, at the same 
time, the control power increases, 


The bearings for the clutches require very careful design with seals to prevent the escape of the mixture, 
especially in the presence of heating. 


The nonlinearity of the clutch's static characteristic and the impossibility of achieving reversal with one 
clutch necessitate the application of special arrangements for switching the clutches in control systems, thus 
providing, in particular, rectification of the characteristic and the production of a reversal. 

SUMMARY 


1. The properties of ferromagnetic clutches make them very suitable for use as actuating devices in fast- 
response automatic control systems, 


2. The starting time and control power for ferromagnetic clutches are considerable smaller than for other 
arrangements (e.g., electric motors). On the basis of a, the ferromagnetic clutches are significantly superior to 
electric motors, friction clutches, and, in some cases, hydraulic actuating devices, The negligible control power 
required for ferromagnetic clutches permits the use of a fast-response control unit (for example, an electronic 
one) having a small time constant, which leads to an increase in the quality of the system. 


3. The response time of ferromagnetic clutches greatly depends on their construction, The reduction of the 
response time requires means of ensuring a rapid build-up of the magnetic flux, as well as means of preventing 
the possible migration of iron particles between gaps, under the influence of centrifugal forces. 


In conclusion the author expresses his gratitude to A. A. Feldbaum for his comments on the paper and his 
help. 
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ADDITIONS TO THE TABLE OF OPTIMAL CHARACTERISTICS 
GIVEN BY SOLODOVNIKOV AND MATVEEV 


N. A. Smyrova 
(Leningrad) 


In a paper by V. V. Solodovnikov and P. S. Matveev [1], devoted to the synthesis of corrective devices for 
servo systems on the basis of given requirements on the dynamic accuracy in the presence of noise, there was 
given a table of optimal characteristics. Expressions were given in this table for the optimal impulse transient, 
k(t), and transfer, @(s), functions, and also the mean-square error for different spectral densities of the random 
components of effective signal and noise, when certain first coefficients of error, Cy, Cy, Cg, ..., are given. The 
gradual accumulation of computed results allows this table of optimal characteristics to be expanded, which 
might be quite useful for calculations. 


This paper presents computed results which supplement the table of the paper cited. 
We shall consider the case when the spectral density of the noise has the form: 


243? (w? +- a? + «?) 
cot + 2 (a? — w®) wo? + (a? + az" (1) 





S, (@) = 


The spectral density of the random component of the effective signal, S,,(w), equals zero, and the dynamic 
accuracy of the reproduced effective signal is determined by given values of the error coefficients, Cy, Cy, C,, 
Cs and Cg, This last condition corresponds to an effective signal in the form of a fourth-degree polynomial in 
time (r = 4). 


In this case, the optimal transient function has the form: 


k(t) = Ag + At + Agt® + Agt® + Agt*+- Bye’ + Bye *T + 
+ K18(t) + Ds(¢—T) for Ot KT, 
k(t)=0 for t>T, (2) 
where } = a* + «2. 


The optimal transfer function is defined by 
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The coefficients Ay, Ay, Ag, Ag, Ag, By, By, Ky and Dy, entering into Formulas (2) and (3), are the solutions 
of the following system of nine equations: 
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The mean-square value of the noise at the output of the optimal system is determined by the expression 
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(13) 


Using the optimal characteristics provided by (2)-(13), corresponding to the case when r = 4, we can 
transfer our attention from them to the cases when r = 3, 2 or 1, i.e., to the cases when the effective input 
signal is expressed by a polynomial of the third, second or first degree in time. 


To do this, it is sufficient to carry out the following: 

1) for r = 3, in Formulas (2)-(13) set Ay = 0, and discard Equation (12) from the system (4) through (12); 

2) for r = 2, in Formulas (2)-(13) set Ay = As = 0, and discard Equations (11) and (12) from the system (4) 
through (12); 

3) for r = 1, in Formulas (2)(13) set Ag = Ag = Ag = 0, and discard Equations (10), (11) and (12) from the 
system (4) through (12). 

It is easily verified that the cases r = 1 and r = 2 give the same results as are to be found in the referenced 
paper [1] (Cf. items VII and VIII). 


Thus, Formulas (2)-(13) supplement the table of paper [1], and provide the possibility of computing 
optimal system characteristics when the effective signal is expressed as a third-degree or fourth-degree polyno- 
mial in time, and the noise . spectral density has the form givenin (1). Moreover, the example considered 
here shows that, for setting up a table of optimal characteristics, it suffices to provide complete data only for the 
effective signal with the most complex form (in the given case, with r = 4), and to state, for the remaining 
cases, which of the coefficients are to be set equal to zero, and which of the equations are to be discarded. 
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